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NOMENCLATURE 


Definition Uni ts 

2 

Flow area between primary nozzle exit and secondary M 

chamber wall as defined in Figure 31 

2 

Secondary chamber throat area M 

Crocco number defined in Section IV.B.3 

Mass flow entrained on the base flow side of the Kg/s 

dividing streamline 

Total mass flow Kg/s 

« 

Length' between primary nozzle exit and secondary M 

chamber throat as defined in Figure 17 

Distance along shear layer of plume M 

Mach number of the plume boundary (Figure 35) 

Mach number at the exit of the primary nozzle 


2 

Static pressure of plume boundary N/cm 

2 

Primary chamber stagnation pressure N/cm 

2 

Secondary chamber pressure N/cm 

Radius of primary nozzle plume used to define plume M 

shape in terms of r/r* 

c 

Radius of primary nozzle exit used to^ define plume M 

shape in terms of X/r^ 

Primary chamber throat radius M 

Static temperature of the plume boundary °K 

Stagnation temperature of primary stream °K 

Axial distance from primary nozzle exit used in M 

defining plume shape 

Dimensional distance measured from the centerline of M 

the shear layer 
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NOMENCLATURE (cont.) 


Term Definition 

' Y Expansion coefficient (ratio of heat capacities - 

6 Plume expansion angle, defined as the nozzle, lip 

angle Gn plus the Prandtl -Meyer expansion angle 
Av in Section IV. B. 2 (Figure 35} 

<5* Displacement thickness across the shear layer as 

defined in Section IV.B.3.C. 

2 

E Primary nozzle area ratio (Rep/Rtp) as described 

P in Figure 3 

2 

e Secondary nozzle area ratio (Res/Rts) as described 

in Figure 3 

6c Momentum thickness across the shear layer as defined 

in Section IV.B.3.C 

a Shape factor of the shear layer 

(j) Velocity profile 


Units 


degree 
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I. 


SUMMARY 


A. STUDY OBJECTIVES AND SCOPE 

The major objectives of this program were to: (1) evaluate the 

results of dual throat thruster cold flow tests, generated on a joint NASA- 
MSFC and ALRC effort, (2) conduct additional tests where required to complete 
the data matrix, (3) define the operating characteristics of the MSFC dual 
throat thruster concept, and (4) generate a computer model that adequately 
represents the dual throat aerodynamics and allows the prediction of per- 
formance for dual throat configurations. 

To accomplish the objectives, a four task program, summarized in 
Figure 1, was conducted. 

Test data were analyzed for a total of 69 tests conducted with 
the dual throat hardware supplied by the Marshall Space Flight Center. The 
first 40 tests were conducted under ALRC IR&D funding. The configurations 
tested on the IR&D program are shown in Figure 2. An additional 29 Mode II 
plume attachment tests were completed on this contract. The Mode II hard- 
ware configurations tested are illustrated in Figure 3. 

An aerodynamic mixing model was formulated to correlate the 
test data. The logic used to develop the aerodynamic model and to incor- 
porate it into a performance model for dual -throat engine performance is 
illustrated in Figure 4. The cold flow test data was used to gain basic 
understanding of dual throat operating characteristics and to Identify the 
most important design concept variables. The data were also used to iden- 
tify the influencing physical mechanisms such as shear mixing, plume 
shape and impingement, weak compression shocks, and secondary wall two- 
dimensional expansion. Identification of the important physical mechanisms 
led to selection of analytical methods of problem solution. 

The methodology was established for predicting the performance of 
dual throat thrusters and for making a preliminary performance prediction 
for a hot firing type test. This methodology incorporated the aerodynamic 
model and utilized existing simplified JANNAF models for performance pre- 
diction. 


A preliminary engine system model was defined utilizing the 
design criteria generated in this program. 

B. RESULTS AND CONCLUSIONS 

The most Important conclusion from this study is that the dual 
throat concept is aerodynamically feasible. The performance loss obtained 
during cold flow was as low as 0.5 percent. The loss is minimized by an 
optimum nozzle spacing corresponding to an Af/Ats ratio of about 1.5, or an 
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_ Figure 1, Dual Throat Thruster Cold Flow Analysis Program Logic 
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Figure 2. Dual Throat Test Configurations 
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Figure 3. MSFC Contractual Test Series, Dual Throat Cold Flow 
Mode II Hardware Configurations 



Figure 4. 


Dual Throat Performance Model Development 








I, B, Results and Conclusions (cont.) 


Le/Rtp ratio of 3 for the 'dual -throat hardware tested. Bleed flow require- 
ment at this condition is only 4 percent. The results indicate that the 
Mode I and Mode II geometry requirements are compatible and pose no signi- 
ficant design problems. 

C. CONTRIBUTION TO THE DEFINITION OF HIGH PERFORMANCE ENGINES 

This study has been successful in confirming, from a performance 
viewpoint, that the dual throat thruster concept can be competitive with 
other tripropellant and mixed mode engines being considered for single- 
stage-to-orbit (SSTO) missions. It remains for these study results to be 
used in an engine system analysis to determine the performance/weight trade- 
offs and the vehicle/engine integration requirements. 
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II. 


INTRODUCTION 


A. BACKGROUND 

Based on both Apollo and Space Shuttle experience, the lead 
time from an "Authority to Proceed" to an operational system is of the 
order of eight to ten years, and the planning for, and the development of 
the necessary technology base for such programs precedes this time by ten 
to fifteen years. Utilizing this knowledge, NASA has made the assumption 
that a follow-on to the present Space Shuttle system will be required in 
the 1995 time frame, and since 1973 has conducted studies directed toward 
the examination of potentially attractive earth-orbit transportation 
systems as a means of identifying the associated technology requirements. 

Advanced high pressure bi propel! ant engine and tri -propel! ant 
qngine (TPE) concepts have emerged from these studies. These propulsion 
systems meet the unique requirement of high performance density demanded by 
the single-stage-to-orbit and heavy lift launch vehicle concepts, and thus 
offer a means to achieve an economical space transportation system. 

One of the major propulsion system candidates that is being 
studied is the dual throat concept, proposed earlier by NASA-MSFC to obtain 
a large area ratio adjustment within a single thrust chamber assembly, 
without the need for extendible nozzles. The concept is readily adaptable 
to TPE applications, and preliminary analysis has indicated a high nozzle 
efficiency. 


B. PURPOSE AND SCOPE 


The feasibility of an economical space transportation system is 
heavily dependent on the delivered performance of the engine system. The 
first step in evaluating a new engine concept for this application, there- 
fore, is to determine the performance. It is the purpose of this study to 
correlate the results of nozzle cold flow tests conducted with a NASA-MSFC 
designed dual -throat thruster at the Aerojet Liquid Rocket Company, and to 
provide a means for determining the performance of a dual-throat engine. 

C. GENERAL REQUIREMENTS 

For purposes of this program both parallel burn and series 
burn conditions were assumed as requirements for the engine during Mode I 
operation. Mode I of the mixed-mode propulsion concept involves the 
sequential or parallel use of high density impulse propellants and high 
specific impulse propellants in a single stage to increase vehicle per- 
formance and reduce vehicle weight. The concept has been described in the 
literature (e.g., see Reference 26)and theoretically provides the per- 
formance necessary for a VTOHL-SSTO Space Transportation Vehicle. 
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II, C, General Requirements (cent.) 


Mode II operation is assumed to involve only the high specific 
impulse propellant. 

D. APPROACH 

To accomplish the program objectives, an effort involving 
four technical tasks was conducted. Tasks accomplished were: 

1. Task I: Dual Throat Cold Flow Test Data 


Additional tests were conducted and the nozzle cold flow 
test data were analyzed to determine the performance potential of the dual- 
throat thruster, 

2. Task II: Aerodynamic Mixing Model 

Significant theories concerning plume shape and plume 
attachment were reviewed and compared with the results of the cold flow 
test data. An aerodynamic mixing model was formulated, a computer program 
was written, and the program was utilized to correlate the test data. 

3. Task III: Performance Methodology 

A method for predicting the performance of dual -throat 
thrusters was established and a preliminary performance prediction for a hot 
firing type test was made. 

/ 

4. Task IV: Preliminary Engine 'System Model 

A preliminary engine system model was defined utilizing 
the design criteria generated on this program. 


8 



m. DUAL-THROAT COLD FLOW TEST DATA 
A. OBJECTIVES AND GUIDELINES 

A total of 69 tests were conducted with the dual throat hardware 
supplied by the Marshall Space Flight Center. The tests were divided into 
three series: (1) plume attachment; (2) secondary only flow, and (3) parallel 

flow. A pictorial description of the tests and corresponding objectives are 
shown in Figure 2. Dual throat terminology is illustrated in Figure 5, which 
also includes a comparison between typical high pressure TPE design parameters 
and those for the dual throat cold flow test hardware. 

The major objectives of this study are to; (1) substantiate 
previous analytical results, (2) provide data concerning dual throat per- 
formance, and (3) correlate the test- results to actual engine design 
conditions. 


B. EXPERIMENTAL APPARATUS 
1 . Dual Throat Assembly 

The dual throat assembly supplied by Marshall Space Flight 
Center consists of two concentric nozzles, spacers for changing their rela- 
tive axial positions, an injector for hot firing the inner nozzle, and mani- 
folding for injecting bleed flow into the annulus between the inner (pri- 
mary) and the outer (secondary) nozzles. The secondary nozzle is equipped 
with 48 pressure taps in 21 axial positions. Thirty-six of the 48 avail- 
able pressure taps were used during this test program. 

The Dual Throat Assembly was modified to allow the high 
volume flows required for cold flow testing with gaseous nitrogen. To 
increase the flow capacity of the primary circuit the injector was removed; 
a new end cap with a two inch diameter inlet pipe was fabricated; and a 
sleeve was installed between the end cap and the chamber to fair the inside 
diameter of the inlet pipe to the inside diameter of the chamber. The capa- 
city of the secondary circuit was increased by removing the diffuser ring, 
drilling a 1.031 cm (.406 in.) diameter hole between each of the .3175 cm 
(.125 in.) diameter holes in the secondary manifold, and replacing the four 
1.588 cm (.625 in.) tubes on the secondary manifold cover with six 2.54 cm 
(one-inch) diameter tubes. Figure 6 shows the configuration of the dual 
throat assembly as it was tested. 

The primary nozzle was modified late in the program to 
allow testing with two primary nozzle area ratios, and to verify the back 
pressure measurement obtained initially with secondary wall pressure measure- 
ments (cf. Figures 7 and 8). Six additional pressure taps were added to the 
primary nozzle outer wall. These measurements agreed well with the secon- 
dary wall measurements, substantiating the early test data. 
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Figure 5. Dual Throat Terminology 




• DIRECT MEASUREMENT (USED IN ALL THREE TEST SERIES) 

• PRESSURE - SECONDARY WALL, PRIMARY SECONDARY P^ 

• TEMPERATURE - PRIMARY AND SECONDARY INLET 

• FLOW RATE - SONIC NOZZLE + REQUIRED PRESSURES & TEMPERATURES 




Figure 6. Dual Throat Assembly as Tested 
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Figure 7. Dual Throat Nozzle Extension 
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Figure 8. Primary Hardware Modifications 







Ill, B, Experimental Apparatus (cent.) 


2. Test Facility 

Figure 9 shows a schematic diagram of the test set up. The 
major components were the gaseous nitrogen storage tank, tank safety valve, 
pressure regulators, pebble bed heat storage device, electric heater, flow 
measurement nozzles, thrust chamber valves, and the test engine. Two flow 
circuits were provided, a low flow circuit with a run line fabricated from 
3.81 cm (1.5 in.) stainless steel tubing, and a high flow circuit which 
used 5.08 cm (two-inch) stainless pipe upstream of the flow measurement 
nozzle and 7.62 cm (three-inch) pipe downstream. The high flow circuit 
contained a pebble bed heat storage device which heated the nitrogen to 
prevent condensation shocks from occurring in the nozzle. The pebble bed 
could maintain the temperature at the inlet to the engine at 355.6 + 18®K 
(180 + 10°F) for 15 sec, at a flow rate of 4.54 Kg/s (10 Ib/sec). Figures 
10 and 11 show front and rear views of the dual throat assembly on the test 
stand. 


3. Instrumentation 


The engine was instrumented for measuring the flow in both 
primary and secondary circuits, gas temperature and pressure in both circuits, 
primary and secondary chamber pressure, and 36 wall pressure measurements 
connected to two scani valves. Each scani valve consisted of a pressure trans- 
ducer close coupled to a solenoid operated pneumatic switch. The pneumatic 
switch had 48 input ports and the system was capable of stepping from one 
port to the next at rates up to 30 steps per second. For this investiga- 
tion two scani valves were used; one for low pressure and one for high pres- 
sure. The pressure taps were connected to the scani valves with plastic 
tubing. Each pressure tap was connected to two scanivalve ports so that 
two separate readings could be obtained each time the scanivalve made a 
complete cycle (see Figure 12). An ambient pressure and a reference pres- 
sure were also connected to each scanivalve. Table I lists the instrumen- 
tation. 


The heart of the instrumentation system was a mini-computer. 
This computer controlled the test, and gathered and processed the test data. 

Two special computer programs were written for the dual throat testing. 

"APLME" was a special data reduction program and "ASCAN" was a modified 
sequencing program. ASCAN performed the normal sequencing functions such 
as operating the thrust chamber valves and turning on the oscillograph at 
operator selected times, checking for malfunctions, and acquiring data and 
converting it to engineering units based on calibration information; but had 
the additional capability of operating the scanivalves. The number of pressure 
measurements connected to each scanivalve, number of data records obtained, 
and the time between completion of a scanivalve step and reading the pressures 
were input variables to the program. 
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Figure 10. Dual Throat Cold Flow Test Stand Front View 





Figure 11. Rear View of Cold Flow Test Stand 
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Figure 12. Secondary Chamber Pressure Instrumentation 


Ill, Dual -Throat Cold Flow Test Data (cont.) 


C. TEST PROCEDURE 

Three test series (40 tests) were conducted under ALRC IR&D 
funding. The configurations tested are shown in Figure 2. An additional 
29 Mode II (Series I A) tests were completed on this contract. The hardware 
configurations for these tests are illustrated in Figure 3. Cold flow 
test history is summarized in Figure 13. 

For Series I and lA (Plume Attachment Tests) the high flow cir- 
cuit was connected to the primary chamber. The primary flow rate was 4.54 
kg/s (10 Ibm/sec) for the entire series, and the bleed flow was varied from 
zero to 0.681 kg/s (1.5 Ibm/sec) for various primary nozzle positions. 

For Test Series II and III the high flow circuit was connected 
to the secondary chamber. A secondary flow of 6.36 kg/s (14 Ibm/sec) was 
used for Test Series II (Secondary Blockage Tests). There was no flow 
in the primary circuit for this series and five nozzle positions were used. 

Test Series III (Parallel Flow Tests) was conducted at a single 
nozzle position with a total flow of 6.81 kg/s (15 Ibm/sec) and primary flow 
fractions of 12 to 29%. For this series both primary and secondary flow 
leads were investigated. 

All test series were conducted according to the following basic 
procedure. First the electric heater was turned on and heated nitrogen was 
allowed to flow through the pebble bed heat storage device. Approximately 
one hour was required to heat the pebble bed. Two seven second tests could be 
run on each heating. While the pebble bed was heating, the primary nozzle 
was adjusted to the proper position as specified in the test plan, and final 
instrumentation calibrations were obtained. When the pebble bed had reached 
367°K (200°F) the electric heater was turned off and the preheat valves 
closed. The tank safety valve was then opened, and the pressure regulators 
set to give the required flow in each circuit. The arm and fire switches 
were then depressed which turned control over to the mini -computer. One 
half second after the fire switch was pushed, the computer turned on the 
oscillograph and the digital data system. Another half second later the 
thrust chamber valves were opened. (For Test Senes III the valves were 
opened sequentially two seconds apart, for Series II only the high flow 
valve was opened, and for Series I the valves were opened simultaneously.) 

The pressures were allowed to stabilize for two seconds after opening the 
valves, before the computer called for the first data record. For each 
data record all important parameters were recorded including one pressure 
measurement for each scanivalve. After each record was obtained the scani- 
valve was advanced one step. When the computer received a signal from the 
scanivalve indicating the step had been completed and the scanivalve was in 
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ALRC IR&D TESTS 


- 40 TESTS 

MODE II - PLUME ATTACHMENT TESTS 
4 NOZZLE POSITIONS 1 
0-10% BLEED FLOWS J~ 

- 21 TESTS 

MODE I - SECONDARY BLOCKAGE 
5 NOZZLE POSITIONS 

- 5 TESTS 

MODE I - PARALLEL BURN 
NASA-MSFC CONTRACTUAL TESTS 


- 14 TESTS 

MODE II - PLUME ATTACHMENT TESTS 

VERIFICATION OF 2 NOZZLE 
POSITIONS 

HARDWARE CONFIGURATION 
VARIATIONS 

2 PRIMARY NOZZLE AREA ' 

- 29 TESTS 

- 9 TESTS 

RATIOS 

2 SECONDARY NOZZLE' 

- - 20 TESTS 

AREA RATIOS 

J 



Figure 13. Dual Throat Cold Flow Test History 



Ill, C, Test Procedure (cont.) 


position' for the next pressure reading, the computer waited 80 milliseconds 
for the pressure reading to stabilize before calling for another data record. 
(Prior to starting the test program, the scani valve was connected to two 
known pressure sources and stepped from one to the other while the output 
signal was recorded on an oscillograph. The results indicated a delay of 
80 milliseconds was sufficient to insure accurate pressure readings.) After 
stepping through the scani valve channels, the computer terminated the test 
by closing the thrust chamber valves. The computer then printed a summary 
of the test data and transferred the average value of each parameter to the 
data reduction program "APLME". APLME first corrected all the scani valve 
measurements based on the reference pressure readings, and tabulated the 
results according to axial distance from the inlet of the secondary chamber 
as shown in Table II. (This was effectively calibrating the scanivalve 
transducers during each test run, and eliminated any errors due to instru- 
mentation drift.) APLME next calculated the flow rate for each circuit and 
for Test Series I also computed the thrust coefficient and specific impulse 
from the pressure measurements. A sample printout is given in Table III. 

D. DATA ANALYSIS 

1. Series I and lA - Plume Attachment Test Results 


The wall pressure integration technique shown in Figure 14 
for calculating thrust is estimated to have an accuracy of + 0.5 percent. 

For this reason and in an effort to reduce program scope, no direct thrust 
measurements were made. As indicated in the figure, the total engine thrust 
is the summation of the primary nozzle thrust (which was obtainedusing the 
thrust coefficient equation) and the integrated wall pressure thrust. The 
engine specific impulse is the total engine thrust divided by the total 
gas flow rate (primary flow plus bleed flow). 

A procedure was developed to isolate the losses associated 
with the dual throat engine concept, i.e., the plume attachment losses and 
bleed flow losses. The test one-dimensional engine specific impulse (IspODE^ 
was calculated as the mass average of the primary and bleed flow circuits 
IspODE values. The IspnpE each circuit was based upon a calculation for 
expanding nitrogen gas to an area ratio of 20 to 1 using the respective 
circuit gas chamber temperatures. (The dual throat hardware Mode II area 
ratio was 20:1.) The percent Tsp loss due to the dual throat concept was 
then taken to be the IspxoK ® conventional nozzle with the same conical 
half angle as the dual throat hardware (0 = 12.5°) minus the measured test 
Ijp divided by the Ispnnr percent Isp loss due to the viscous 

drag within the boundary 1 ayer of the nozzle at the cold flow test condi- 
tions. This loss was calculated using the TBL chart program and was found 
to be less than 0.1% and, therefore, was ignored in this analysis. This cal- 
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TABLE II 


CORRECTED WALL PRESSURES 

dual throat ffAftOWARf ^2 COCB FLOI'I TF.ST program 

TEST AJO. CFR6-7P7-I09 

TEST OATE 02-0«-7'8' R.SO.UJ' 

BAROMETRIC PRESSURE 1U.A63 

REFERENCE PRESSURE .01930 .01930 “ 


SECTION 3 — secondary CHAMBER WALL PRESSURE FROM SCANIVALVE INPUT 


HO 

CO 


X 

R 

TIME 

PWl 

PCH 

PWl/PCH 

TIME 

PW2 

PCH 

PW2/PCH 

TIME 

pv*3 

PCH 

PW3/PCH 

CIN) 

CIM 

(9FC) 

(PSIA) 

<PSIA) 

.. , 

(SEC) 

(PSIAJ 

(PSIA) 


(SEC) 

(PSIA) 

(PSIA) 

_ 

2.900 

2.925 

2.089 

11.650 

323.815 

.0360 

3.150 

1 1.669 

320.806 

.0359 

.000 

,000 

.000 

. 0000 

0.037 

2. 672 ' 

2.239 

11.700 

320.053 

.0361 

.000 

.000 

.000 

.0000 

,000 

.000 

.000 

.066TI 

5.190 

2.621 

2.361 

11.726 

320.370 

.0361 

3.300 

11.759 

325.375 

.0361 

.000 

.000 

.000 

,0000 


5.939 
6.029 
6. '90 0 
7.169 
7.000 


2.197 
1.971 
1.637 
1 .806 
1.795 


2.097 

2,069 

2.239 

2.361 

2.097 


11.709 
1 1.596 
11.013 
11.305 
11.301 


320.688 
323.815 
320,053 
320.370 
320. 688' 


.0361 

.0358 

.0352 

.0309 

.0308 


3.022 
.000 
3.681 
_ .000 
3.'630 


11.750 

.000 

11.519 

_.000 

11.372 


325.719 

326.033 
,_p00 
T?6 .oo'6 ■ 


.0000 
'.0353 
. 0000 
.030 8 
. 0000 


.000 

■flT36l 

.000 


.000 


11.007 

.000 


0.083 

.000 


10.903 

.000 


_ _i,ooo 

32'7.7Sa" 

. 000 

328.098 

.000 


.0000 
To 308 
.0000 
.0332 
.0000 
.0268 
.00 00 
“.ooo^Oq 

. 0 liT 
.0000 


7.940 

1.637 

2.769 

8.970 

324.529 

.0276 

3.953 

8.790 

326.406 

8.4A1 

1 .644 

2.692 

6.272 

324.820 

.0163 

4.089 

6.263 

327.120 

6.940 

2.055 

3.028 

5.073 

324.688 

.0156 

.000 

,000 

.000 

9,440 

2.166 

3.150 

4.380 

324.846 

.0135 

.000 

.000 

.000 

9.939 

2.276 

3.300 

3.998 

325.375 

.0123 

4.2! 1 

3.969 

327.437 

10.436 

2.367 

3.422 

3.365 

325.719 

.0103 

.000 

.000 

.000 


.0269 
.0 191 

.onoo 

.0000 

.0121 


0.619 

,000 


.000 

.000 


6.791 

_.00J) 

. 000 “ 

.000 


0,701 


3.832 
.000 


328.072 
._0(^ 
.(iOO 
.000 
326.508" 
.000 


10.930 

11.437 

11.937 

12.436 

12.927 

13.508 

2.497 

2.608 

’2.719 

2.630 

2.939 

3.067 

3.558 

2.620 

2.769 

.000 

2.892 

3.028 

2.836 

2.466 

2.075 

.000 

1.604 

1.315 

326.300 

324.688 

324.52"9 

.000 

324.820 

324.680 

.0087 

.0076 

.0064 

.0000 

.0049 

.0041 

.000 

,000 

3.558 

3.661 

.000' 

.000 

,000 

.000 

2.i42' 

1.781 

.000 

.000 

.000 

.000 

‘326.300’ 

326.433 

,000 

.000 

.0000 
.0000 
.0 066 
.0055 
'. oo’oo 
.0000 

.000 
.000 
3.953‘ 
.000 
’ ^’.’o’oo 
.000 

.000 
.000 
2’.106’’ 
.000 
.000 " 
.000 

.000 

.000 

'326. 406" 
.000 
■ .000 
.000 

.0000 
.0000 V 
■’.’0665’ , 

.0000 

.0000 

.0000 

2.940 

2ND SCAN 
2.925 

4.077 

11.776 

328.785 

.0356 

5.936 

11.820 

329.7’64 

.0358 

• O'OO 

.000 

.000 

.0000 

4.437 

2.672 

5.013 

11.630 

328.891 

.0360 

.000 

.000 

.000 

.0000 

.000 

.000 

.000 

.0000 

5.190 

2.621 

5.149 

11.637 

329.024 

.0360 

6.061 

11.902 

329.975 

.0361 

.Tfto’o 

,060 

.000 

.0000 

5.939 

2.197 

5.27? 

11.633 

329.261 

.0359 

6.197 

11.887 

330,689 

.0359 

6.591 

1 1,886 

331.376 

.0359 

”6.42'9 

’1.971“ 

4.8T7 

1 1 :7TS“' 

328.785 

“ .“OTSTS" 

.000 

,000 

TWO 

.0006 

.000 

.000- 

.006 

.0000 

6.940 

1.837 

5,013 

11.551 

326.891 

.0351 

6.469 

1 1 .656 

330,557 

.0353 

7.121 

11.535 

331.800 

.0348 

7.189 

1 .806 

5.149 

1 1.44? 

329.024 

.0340 

.000 

.000 

.000 

.0000 

.660 

.O'OO 

.000 

,5000 

7.440 

1 .795 

5.272 

11.350 

329.261 

.0345 

6.591 

11.580 

33! .376 

.0349 

7.257 

11,050 

331.905 

• u ^ J i 

7.691 

1.006 

5.408 

10.627 

3?9.286 

.0323 

’. 000’ 

.’OOO 

■ ■■ . O'OTi 

."5oW 

.000 

.000 

.000 

.Trno6 

7.940 

1 .637 

S.S30 

9.122 

329.420 1 
— 7' •50 — •rrm — 

.0277 

A ro i — 

6.727 

L 'a/ia — 

8.968 

X Cti“' 

331 .958 

.0270 
TIS 10? 

7.380 
. nnO 

8.933 

rtrtrc 

332.355 
nnnr" 

.0269 

— nnnnr" 



6.900 

9,000 


2.055 

2.166 


5.002 

5.938 


5.156 

0.053 


329.073 

3?9.760 


,0156 

.0135 


.000 

.000 

6.985 


.000 

.000 

0.016 


.000 
.000 
331.376 


.0000 

■.WOTT' 

.0121 


'.T5W 

7.516 


■,froo” 

3.669 


" .THTO- 
332.619 


.0000 

,0116 


10.438 

10.938 

2.387 

2.497 

6.1R7 

6.319 

3.411 

2.870 

330.689 

330.716 

V0 103 
.0087 

‘.000 

,nno 

“.OOO" 

.000 

■ TOD 0 
.000 

■ loooO 
.0000 

•."000“ 

.000 

TCOTJ" 

.000 

.“OTTO' ■ 
.000 

■.06 06 

.0000 
■ A A A A 

ri'.437“ 

2.608 

5TW8 

-g-.'Wir’ 

329;2'8B — 

.0076 

TOTIO 

nroiT 

.00 0 

.0000 

,000 

,000 

.000 

• Ofi 0 0 

1 1.937 

2.719 

5.530 

2.103 

329.420 

.0064 

6.319 

2.166 

330.716 

.0066 

6.727 

2.147 

331.956 

.0065 

12.416 

2.830 

.oOO 

.000 

.000 

. 0000 

6.469 

1 .807 

330.557 

.0 055 

.00D‘ 

.000 

.000 

.0600 

12.927 

2.939 

5.680 

1.604 

329.790 

.0049 

.000 

.000 

.000 

.0000 

.000 

,000 

.000 

.0000 

13.508 

3.067 

5.80? 

1.291 

329.473 

.0039 

.000 

.000 

.‘OO’O * 

.“0060" 

.065 

~ ■:(roo 

■ " .TOO 

“ .’O'OOO 



TABLE III 


FLOW AND PERFORMANCE CALCULATION (1 of 3) 


DUAL Throat gn2 

TEST NO, CFBS 
TEST date 02 

TYRE w 1 

COLO PLOW TEST 
•727«i09 , 
-08.78 12,50,43 
program input 

PROGRAM 


type 

type 

1 9 MODE 2 

2 P MODE 1 

PLUME attachment 
SERIES AURN test 

TEST 


type 
uertp p a 

3 ^ MODE 1 PARALLEL BURN TEST 
(.883 

primary plow RATE (LBM/8EC 

10,004 


SECONDARY FLOW RATE (LBM/SEC) 

,390 


primary static chamber PRESSURE (PSIA) 
PRIMARY GAS temperature CR) 

328,406 

663,629 


secondary 6A8 TEMPERATURE (R) 
average of wall Pressures (psia) 

SU,73S 

p 

X (INCHES) 

R (INCHES) 

PRESSURE 


1 

a'9A0 


11.706 


2 

«.«37 

2.872 - 

11,855 


3 

5j.l90 

2,621 

11,850 


M 

B.959 

2,197 

11.815 


5 

6,429 

1,971 

11,731 


6 

6.940 

1,857 

11,519 


7 

7.189 

1,006 

11.435 


a 

7,440 

1,795 

11, 2S? 


9 

7,691 

1,806 

10,574 


10 ' 

7.940 , 

1,837 

8,919 


u 

8,441 ' 

1,9«4 

6,515 


12 

8,940 

2,055 

5,155 


15 

9,440 

2,166 

4.431 


1« 

9.939 

2,276 

3,951 


15 

10,438 

2,3S7 

3,390 


16 

10,958 

2,497 

2,852 


17 

11,437 

2,608 

2,490 


18 

11,937 

21719 

2,125 


19 

12,436 

2,830 

1.794 


20 

12,927 

2,939 

1,609 


21 

13,508 

3.067 

1,309 



average primary NOZZIE OUTER WAlU PRESSURE 

PnBaR • 11.81T PSIA 






dual throat GN2 COLD FLOW TEST PROGRAM 
TEST Nn. CFBfe-727-109 
TEST DATE 02-06-78 i2,30,a3 
PROGRAM INPUT 

TYPE R 1 

type I m MODE 2 PLUME ATTACHMENT TEST 
TYPE 2 » mode i SERIES BURN TEST 
type 3 ■ mode 1 parallel burn TEST 
LERTP Q 2.685 

PRESSURE calculations 

STATIC STAGNATION 
PRIMARY jSa'^406 339*520 

SECONDARY 11.786 11*809 

POS/POP a .035 
FLOW CALCULATIONS 

PRIMARY SECONDARY 

FLOW RATE (LB/SEC) lO.OOii ',590 

MAX’, IDEAL flow lO.lttS 2;?67 

percent of max, 98,598 21,325 

normalized flow ,759 1,131 

Flow coefficient ,9S6 ,2i3 

X bleed flow 9 5',570 
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TABLE III 


FLOW AND PERFORMANCE CALCULATION (3 of 3) 


DUAL THROAT 6H2 COLO ^LOM TEST PROOPAM 

Test no. cpbs-tzt-io# 
test date 02-06-78 12. JO, <13 
PR09RAH INPUT 

TYPE • J 

type i m HOOE 2 PLUME ATTACHMENT TEST 
TYPE 2 « MODE 1 SERIES BURN TEST 
type J ■ MODE I parallel BURN TEST 
LERTP 9 2,083 

THRUST CALCULATIONS 

PRIMARY EXPANSION 

primary nozzle thrust (LbP) « 702.732 


nozzle wall pressure XNTEORaTION , 

♦•THRUST generated ALONG PATH FROM PRIMARY NOZZLE EXIT TO SECONDARY CHAMBER, 

normal projected area 9 28,828 IN SQ 

back pressure 9 n,82R PSIA 

THRUST P 288,805 LBPS 

••THRUST generated ALONG RATH FROM SECONDARY CHAH9£R TO SECONDARY NOZZLE EXIT, 


X 

R 

8 

thet 

PM 

DF8E 

DFSE8 

IIN) 

(IN) 

(IN2) 

<pEb) 

tP?lA) 

(L8F) 

<LBF) 

2.980 

21925 

•1965 

>000 

11,786 

•111800 

277,805 

8,837 

2,872 

•8,331 

■10,800 

11.835 

•51,295 

220,110 

5,190 

2,621 

•6,818 

•26,000 

11,850 

•75,938 

150,172 

5,939 

2,197 

•2,939 

•30,000 

11,815 

•38,880 

115.332 

6,829 

1,971 

•1 ,603 

•19 700 

11,731 

•18,636 

96,697 

6,980 

1,837 

-,35s 

•9,000 

11,519 

•8,072 

93,625 

7,189 

1,806 

»,12« 

• 8 800 

11,933 

®!,911 

91,218 

7,880 

1,795 

;i28 

,000 

11,252 

1,358 

92,572 

7; 691 

1,806 

,355 

8,800 

10,578 

5,858 

96,030 

7,980 

1,837 

1,271 

9 600 

8,919 

9,681 

105,711 

8,881 

1,988 

1,595 

12,500 

6,315 

7,988 

113,698 

8,980 

21085 

1,«Z2 

12,500 

5,135 

7,080 

120,735 

9,880 

a, 106 

1,535 

12,500 

9,831 

6,933 

127,168 

9.939 

21276 

1,626 

12,500 

3,981 

5,966 

133,136 

10,830 

213ST 

1,688 

12,500 

31390 

5,268 

138,908 

10,930 

2,897 

1,780 

12,300 

2,852 

8,735 

183,159 

11,837 

2,608 

1,858 

12,500 

2>9o 

8,206 

187,885 

11,937 

2,719 

1,935 

12,500 

2,125 

3,792 

151,237 

12.830 

2>830 

1,975 

12,500 

1,798 

3,361 

158,598 

12,927 

2^39 

2,915 

12,500 

1,009 

3,528 

156,122 

13,506 

3.O0Z 

>38 

121500 

1,509 

1.051 

159,173 


performance summary 


Thrust isp cf ode % loss 

METHOD (LBF) <8ECj (SBC) JSP 

1 801,905 Sl',598 ' 1«72S S2.73T .909 

total BAS PLOW RATB B 10,598 L8M/8IC 
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© PRESSURE INTEGRATION TECHNIQUE, THRUST MEASUREMENT ACCURACY = ±0.5% 



Figure 14. Dual Throat Cold Flow Thrust Calculation 



Ill, D, Data Analysis (cont.) 


culation, as described in Figure 15, normalizes the loss to make it inde- 
pendent of the nozzle contour. 

Plume attachment test data for Series I (IR&D Program) are 
summarized in Table IV. 

The wall pressure profiles for each test, plotted as a 
function of axial distance from the injector face for various amounts of 
bleed flows, are shown in Figures 16 to 19. Superimposed on each pres- 
sure plot is a profile of the two nozzles showing their relative positions. 
For the cases of zero or near zero bleed flows the presence of the plume 
attachment shock is evident by the sudden increase in the wall pressure.. 

As bleed flow is added and the back pressure increases, the plume is 
turned and the attachment deflection angle is reduced to the point where 
the shock is eliminated. 

The corresponding performance data are also shown plotted 
with each pressurepprofile. These dual throat concept performance losses 
are plotted as a function of bleed flow percentage. In each plot, the 
performance losses reach a minimum for bleed flows of 4-8 percent. The 
percent bleed flow corresponding to the minimum performance loss also corres- 
ponds to the minimum bleed flow necessary to eliminate the wall pressure 
spike (plume attachment shock). This trend holds for all four positions 
tested during the IR&D program. A summary of the performance losses of 
the four positions is given in Figure 20. A cross plot of these data is 
shown in Figure 21. These plots indicate that the lowest loss of about 0.5^ 
occurs at an Le/Rtp of about 3.0, with approximately 4.0 percent bleed flow. 
The rapid increase in the performance loss for zero bleed flow for the 
larger nozzle separation distances is a result of the plume attaching up- 
stream of the secondary nozzle throat. As this occurs the plume attachment 
deflection angle is greater and the resultant shock losses are increased. 

The optimum nozzle spacing seems to be obtained when the plume is nmade to 
attach tangent to the wall at the secondary throat plane while using the 
smallest amount of bleed flow. 

A summary of the contractual testing (Series I A) is given 
in Table V. The results of the first 9 verification tests are plotted on 
Figures 18-20 and compared to the IR&D (Series I) results. Agreement is 
good between the two data sets. The apparent bias of 0.5 to I.OiS is 
believed tq be an instrumentation bias as the pressure transducer used to 
measure the upstream pressure of the metering venturi in the high flow 
circuit was ‘decreased from 20.68 x 106 N/m^ (3000 psi) to 13,78 x 106 
N/m2 (2000 psi). The resultant flow measurement uncertainty remained the 
same at + 1%. 
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9 Dual Throat Concept Loss Calculation 
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where: — = 99.3 for a 12.5° conical nozzle 
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Figure 15. Dual Throat Cold Flow Performance Evaluation 


*Percent Isp loss due to boundary layer viscous drag for the conventional 12.5 
degreeVcpnicaT' nozzle 1 (loss .was 0 to. 0.1% for |t|)e^cold flow .test^')- 
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TABLE IV. - DUAL THROAT COLD FLOW SERIES I TEST SUMMARY 
ALRC IR&D MODE II TESTS 


CHAMBER STAGNATION 


TEST 

TEMPERATURE 

°K 

PRESSURE. 

N/ cin^ 

% 

BLEED 

L /R 
e tp 

F 

(KN) 

I^pTEST 

SEC. 

Ve8!' 

LOSS 

PRIMARY 

SECOND 

PRIMARY 

bbUUND 

101 

361 

302 

195 

1 7 

0.0 

1 067 

3 02 

80 7 

82 5 

1 39 

102 

354 

282 

196 

3 3 

4.8 

1.067 

3.17 

79.8 

81.2 

1.04 

103 

366 

286 

203 

5.4 

7.8 

1 067 

3.38 

80 1 

82 2 

1 90 

104 

361 

277 

210 

9.8 

14 2 

1 067 

3.64 

77.9 

80.9 

3.10 

106 

356 

303 

203 

3 6 

0 0 

2 883 

3 13 

79.8 

81.9 

1.90 

108 

363 

291 

204 

6 6 

4.2 

2.883 

3.31 

81 0 

82 3 

0 82 

109, 

356 

286 

208 

8 6 

8 2 

2.883 

3 48 

79 7 

81 2 

1 16 

110 

355 

298 

201 

5 1 

2 0 

2 883 

3 20 

80.4 

81.6 

0 76 

111 

351 

299 

207 

4.8 

1.4 

2 883 

3.26 

79 9 

81 2 

.90 

112 

360 

307 

199 

4.7 

0.0 

5 046 

2 89 

75 9 

82 3 

7.12 

113 

358 

297 

203 

6 7 

4.0 

5.046 

3.19 

78 8 

81 8 

3 01 

114 

362 

288 

203 

8.1 

7 2 

5.045 

3.31 

79.3 

82.0 

2 57 

115 

359 

282 

205 

9.4 

10 3 

5.045 

3.43 

78 1 

81 3 

3.20 

116 

361 

294 

203 

3 9 

3.1 

1 067 

3.22 

79 9 

82 2 

2 15 

117 

356 

324 

206 

3 1 

1 3 

1.067 

3.23 

' 80.1 

81.8 

1 35 

118 

364 

299 

207 

4.0 

2.4 

1.067 

3.28 

80 9 

82 6 

1.38 

119 

358 

304 

209 

1 9 

0 0 

1.067 

3 22 

80.2 

82 1 

1 67 

120 

360 

300 < 

203 

1.2 

0.0 

0.0 

3.12 

79 8 

82.4 

2.46 

121 

355 

285 

206 

3.7 

3.5 

0 0 

3.26 

79.2 

81 5 

2 10 

122 

362 

279 

199 

5.0 

7 5 

0.0 

3.22 

78.4 

31.8 

3.44 

123 

357 

281 

205 

4 5 

5 7 

0 0 

3 29 

78 7 

81 5 

2 74 



TABLE IV (cont.) 



TEWPERATURE. ’R 

CHAK3ER STASNATIOH 
PRESSURE. PS!A 

* 


^ac 

isp test 

'^ODE 

-SEC 

82.5 

xisp 

TEST 

WlflARY 

SECOND 

PRIMARY 

SECOKQ. 

BLEED 

LBF 

SEC 

LOSS 

101 

650 

543 

283 

2 5 

0.0 

1.067 

679 

89.7 

1.39 

102 

637 

508 

CO 

CM 

4.3 

4 8 

1,067 

713 

79.8 

81 2 

1.04 

103 

658 

514 

295 

7 S 

7.8 

1.067 

759 

80.1 

82.2 

1.90 

104 

649 

458 

304 

14 2 

14 2 

1.067 

819 

77.9 

80.9 

3 10 

105 

641 

545 

295 

5 2 

0.0 

2.083 

703 

79 8 

81.9 

1.90 

108 

653 

523 

296 

9.5 

4 2 

2.883 

743 

81 0 

82.3 

0.82 

109 

641 

515 

301 

12.5 

6.2 

2 883 

782 

79 7 

81.2 

1.16 

no 

639 

537 ' 

292 

7 4 

2 0 

2.883 

720 

89.4 

81.6 

0 76 

111 

631 

538 

300 

6.9 

1.4 

2 833 

734 

79.9 

81,2 

.90 

112 

648 

552 

289 

6 a 

0.0 

5.045 

650 

75,9 

82 3 

7.12 

113 

644 

535 

294 

9.7 

4.0 

5 046 

718 

78 8 

81.8 

3 01 

114 

652 

519 

294 

n 8 

7 2 

5.045 

744 

79 3 

82 0 

2 57 

115 

647 

507 

298 

13.7 

10:3 

5.045 

771 

78,1 

81.3 

3.20 

116 

650 

•530 

295 

S 7 

3.1 

1.067 

724 

79.9 

82 2 

2.15 

117 

641 

534 

299 

4 5 

1.3 

1 067 

726 

80.1 

81.8 

1 35 

118 

655 

538 

300 

5.8 

2.4 

1.067 

737 

80 9 

82 6 

1.38 

119 

644 

548 

303 

2 8 

0.0 

1.067 

724 

80.2 

82.1 

1.67 

120 

648 

540 

295 

1.7 

0 0 

0 0 

701 

79.8 

82.4 

2 46 

121 

639 

5i3 

299 

5.3 

3 5 

0.0 

734 

79.2 

81.5 

2.10 

122 

651 

502 

288 

7 3 

7 5 

0.0 

725 

78.4 

81.8 

3.44 

123 

643 

506 

297 

6.5 

5.7 

0 0 

739 

78.7 

81 5 

2.74 
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Figure 16. Dual Throat Cold Flow Mode II Plume Attachment 
Test (Le/Rjp =0.0) 
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Figure 18. Dual Throat Cold How Mode rr Plunie Attachment 
Test (Lg/R^p = 2.883) 
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Figure 21. Dual Throat Cold Flow Mode II Perfo 
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TABLE V. - DUAL THROAT COLD FLOW SERIES I A TEST SUMMARY 
CONTRACT MODE II TESTS 


CHAMBER STAGNATION 


TEST 

TEMPERATURE 

t^RIMARY 

SECOND 

PRESSURE. 

TRTMAR7 

N/cm^ 

SECONir 

% 

BLEED 


c 

*^Vac 

KN 

TEST 

sp 

SEC. 

®PODE 

SEC 

*^sp 

LOSS 

1 


C 

s 

101A 

363 

280 

199 

8.0 

IJ 

5.045 

3 24 

78 5 

81.9 

3.5 

1 

< 1 666 

2 979 

102 

356 

312 

202 

4.6 

0.0 

5 045 

2.90 

74 5 

81.9 

8 2 

\ 




103 

369 

287 

214 

5.9 

2 0 

5 045 

3 22 

77 6 

83 2 

6.0 

\ 




104 

358 

285 

217 

7.0 

4.0 

5 045 

3.39 

77 6 

81.8 

4.4 

1 




105 

370 

279 

223 

8 0 

5,8 

5.045 

3 55 

79 2 

82.8 

3 7 

1 




107 

366 

302 

255 

3.9 

C 0 

2.883 

3.44 

80 2 

83 1 

2.7 

1 




108 

357 

286 

230 

5 7 

2 0 

2 883 

3 63 

80.5 

81 3 

0 8 

i 




109 

368 

284 

234 

8 1 

5 6 

2 833 

3.83 

81 4 

82.7 

1.0 





110 

356 

286 

238 

7 2 

3 7 

2.883 

3 84 

80.4 

81 5 

0.7 





111 

357 

278 

228 

4.8 

8.6 

4 130 

3 70 

77.1 

81 1 

4.3 

1.666 

1.626 

112 

349 

292 

230 

2 2 

0.0 

4 130 

3.41 

76 0 

81 1 

5 6 

S 




113 

358 

282 

209 

3.9 

6.6 

4 130 

3 35 

77.9 

81 5 

3 7 





114 

349 

282 

210 

3 7 

5 7 

4.130 

3 35 

77.4 

80.7 

3.4 





115 

367 

284 

205 

3 0 

3 5 

4.130 

3 21 

79.3 

82 7 

3 5 





116 

365 

285 

208 

2 7 

2.0 

4 130 

3 20 

78.9 

82.7 

3 9 

- 




117 

366 

279 

210 

3 0 ' 

6.2 

3 934 

3 39 

79.5 

82 3 

2.9 

3'. 164 

1.626 

118 

354 

288 

212 

1 1 

0.0 

3 934 

3 25 

79 0 

81.7 

2 6 

i 


t 

119 

369 

279 

208 

2 6 

4 5 

3 934 

3 32 

80 5 

82 9 

2 2 

\ 




120 

363 

281 

211 

2 1 

2.3 

3.934 

3 32 

80.6 

82 5 

1 6 

\ 




121 

368 

280 

196 

2.7 

6 2 

3.213 

3 16 

79.0 

82 6 

3 6 

1 




122 

359 

286 

199 

0.8 

0.0 

3.213 

3.05 

79 2 

82.2 

3 0 





123 

374 

283 

196 

2 3 

4 5 

3 213 

3.12 

80 2 

83.4 

3 2 ' 





124 

366 

285 

200 

1 7 

2.1 

3 213 

3.13 

79 9 

82.8 

2 8 





125 

374 

279 

195 

6.8 

8 0 

4.930 

3.20 

79 8 

82.7 

2 7 

3 1164 

2.979 

126 

364 

285 

198 

4 7 

2 3 

4.930 

3 14 

81.1 

32 6 

1.1 

1 




127 

374 

282 

197 

6.3 

6 1 

4 930 

3 21 

81 0 

83 2 

2 0 

1 




128 

366 

283 

201 

5 7 

4 2 

4 930 

3.23 

81 0 

82.6 

1 2 

1 




129 

376 

277 

198 

7 6 

10 2 

4 930 

3 29 

79 5 

82.9 

3 5 





130 

371 

311 

201 

3 2 

0.0 

4 930 

3.02 

79 9 

83.5 

3 6 

1 






TABLE V. - DUAL THROAT COLD FLOW SERIES I A TEST SUMMARY 
CONTRACT MODE II TESTS (cont.) 


TEST 

TEMPERATURE. °R 
PllltlARY SEtOKD.' 

CHAHOER STAGhUTION, 
PRESSURE, PSIA 
PRIMARY SECOND. 

BLEED 


*^Vac 

LBF 

Isp TEST 
SEC 

^^'’ODE 

SEC 

*isp 

LOSS 





IOTA 

654 

504 

288 

11 6 

7.7 

5 045 

728 

78 5 

81.9 

3.5 

' 1.666 

2.979 

102 

640 

561 

293 

6 7 

0.0 

5 045 

653 

74.5 

81 9 

8 2 





103 

664 

517 

310 

8 6 

2.0 

5.045 

724 

77 6 

83.2 

6 0 





104 

645 

513 

315 

10 2 

4.0 

5.045 

761 

77 6 

81.8 

4.4 





10S 

665 

503 

323 

11.6 

5 8 

5.045 

799 

79.2 

82.8 

3 7 





107 

659 

543 

327 

5.6 

0.0 

2 833 

774 

80.2 

83,1 

2 7 





103 

642 

514 

333 

8.3 

2 0 

2 383 

817 

30 5 

81.8 

0 8 





109 

663 

511 

340 

11 8 

5 $ 

2 883 

862 

81.4 

82.7 

1 0 

t 

f 



no 

640 

514 

345 

10.4 

3 7 

2 883 

863 

80.4 

81.5 

0.7 

\ 




111 

642 

501 

330 

6.9 

S 6 

4 130 

832 

77 1 

81 1 

4 3 

1 ees’ 

1.626 

112 

623 

525 

334 

3 2 

0.0 

4.130 

766 

76.0 

81.1 

5.6 





113 

644 

507 

303 

5 7 

6 6 

4.130 

754 

77 9 

81.5 

3.7 





114 

629 

508 

304 

5 4 

5.7 

4.130 

753 

77 4 

80 7 

3 4 





115 

660 

511 

298 

4 4 

3.5 

4 130 

722 

79 3 

82 7 

3.5 

\ 

f 



116 

657 

513 

302 

3 9 

2.0 

4.130 

720 

.78 9 

82.7 

3.9 

11 


1 

< 

117 

653 

502 

304 

4.4 

5.2 

3 934 

761 

79 5 

82.3 

2.8 

: 3.164 

1.626 

113 

633 

518 

308 

1.6 

0 0 

3 934 

730 

79 0 

81.7 

2.6 





119 

665 

503 

301 

3.6 

4.5 

3.934 

747 

80 5 

82 9 

2.2 

• 




120 

654 

505 

306 

3.0 

? 3 

3.934 

747 

80.6 

82.5 

1.6 





121 

662 

504 

284. 

3 9 

5.2 

3.213 

710 

79.0 

82.6 

3.8 





122 

646 

514 

289 

1.2 

0.0 

3.213 

685 

79.2 

82.2 

3 0 





123 

674 

509 

284 

3.3 

4 5 

3.213 

701 

80 2 

83.4 

3 2 





124 

658 

513 

290 

2.5 

2.1 

3.213 

704 

79 9 

82.8 

■ 2,8 




' 

125 

667 

503 

283 

9.9 

B.O 

4.930 

720 

79.8 

82.7 

2.7 

3.164 

2.979 

126 

655 

513 

287 

6 8 

2.3 

4.930 

706 

81.1 

82.6 

1.1 





127 

673 

507 

286 

9.1 

I5.I 

4 930 

722 

81.0 

83.2 

2.0 





128 

658 

510 

291 

8 3 

4 2 

4.930 

726 

81.0 

82.6 

1.2 





129 

677 

499 

287 

11.0 

'10 2 

4.930 

740 

79.5 

82.9 

3.5 


1 



130 

667 

560 

291 

4.7 

0.0 

4.930 

680 

79.9 

83.5 

3.6 

1 


T 



CO 



Ill, D, Data Analysis (cont.) 


The remaining 20 tests were conducted using different primary 
and secondary nozzle area ratio combinations. The primary objective of these 
tests was to determine the Mode II performance effects produced by different 
geometry configurations and bleed flows. The wall pressure profiles and 
corresponding performance losses are plotted in Figures 22 to 25. 

The same trends in performance loss and shock elimination 
are evident with these data as were evident in the Series I results. The 
summary of these data is shown in Figure 26. These test results show that 
the performance losses are minimized by increasing the primary nozzle exit 
area and decreasing the secondary throat area, indicating that the plume 
attachment shock strength and the resultant performance loss are reduced as 
the gap between the nozzles is eliminated. In the limit, of course, the 
dual throat losses go to zero as the gap is closed forming a conventional 
nozzle. 


2. Series II - Secondary Blockage Test Results 

The performance loss determined by the Series I testing is 
restricted to Mode II operation of the dual throat engine. In practice, 
the compatibility of Mode I and Mode II operation must be considered when 
designing an engine. Therefore, Mode I secondary blockage (Series II) 
tests were conducted during the IR&D effort. 

For all Series II tests 100% of the flow emanated in the_ 
secondary nozzle. The objectives and approach for this series are given in 
Figures 2 and 27. 


Five tests were completed, one with the primary nozzle 
removed, and four other tests with the primary nozzle at discrete positions. 
Table VI is a summary of the tests conducted during this series. 

Figures 28-30 show the measured wall pressure distribution 
along the secondary wall for three primary nozzle positions, respectively. 
Also included in each figure, for comparison, is the wall pressure distribu- 
tion obtained when the primary nozzle was completely removed. 

The tests were run with a constant secondary flow rate of 
about 6.4 Kg/s (14 Ibm/sec). As the primary nozzle is located further from 
the secondary throat, the actual flow area increases thus reducing the plenum 
chamber pressure. 


The secondary blockage test results indicate that blockage 
of the primary nozzle becomes a problem if the separation distance is too 
small and the primary nozzle begins to obstruct the flow in the secondary 
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Figure 22. Dual Throat Cold Flow Mode II Plume Attachment 
Test (L/Rtb = 4.130) 
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Figure 24. Dual Throat Cold Flow Mode II Plume Attachment 
Test (L^/Ryn = 3.934} 
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Figure 27. Series II Secondary Blockage Test (No Primary Flow) Program 



TABLE VI 


SECONDARY BLOCKAGE SUMMARY 



TEST NO. 

SECONDARY FLOW 

^CS 

w 


Kqls 

( Ibm/sec) 

N/CM^ 

(Psia) 


124 

6.21 

03.7) 

25.3 

(36.7) 

Removed 

125 

6.26 

(13.8) 

29.1 

(42.2) 

5.045 

126 

6.21 

(13.7) 

28.6 

(41;54)* 

2.883 

127 

6.40 

04.1) 

43.8 

(63.5) 

1.067 

128 

6.26 

(13.8) 

54.4 ?b.C78.9) 

0,000 


*DATA QUESTIONABLE (DIFFUSER FLOW ATTACHMENT?) 
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Ill, D, Data Analysis (cont.) 


nozzle. To maintain choked flow in the secondary circuit solely at the 
secondary throat plane, the primary nozzle must be located such that the 
flow area between the primary and secondary nozzles is greater than the 
secondary throat area. Figure 31 is a plot of the flow area (Ap) divided 
by the secondary throat area (Ats) versus the nozzle separation distance 
Le/Rtp for both the low (e = 1.666:1) and high (e = 3.164:1) area ratio 
primary nozzles. To maintain a flow area ratio Af/Ats ^ the Le/Rtp 
must be >1.9 for the low area ratio nozzle and >2.5 for the high area 
ratio primary nozzle. The Mode I and Mode II aerodynamic operating 
regions are shown plotted together in Figure 32. 

3. Series III - Parallel Flow Test Results 


In this series of tests both nozzles were operated with 
sonic flow through their respective throats. The flow split (% primary/ 
% secondary) ranged from 10%/90% to 30%/70%. The test objectives and 
approach are given in Figure 33. 

A summary of the parallel tests is listed in Table VII. 
The primary chamber pressure ranged from 34-79 N/cm^ (50-115 psia). The 
secondary chamber pressure from 6.7-41 N/cm^ (10-60 psia). The chamber 
pressure ratio (Pcs/Pcp) ranged from .311 to .912. The primary flow 
ranged from 10-30% of the total flow; the secondary flow ranged from 
70-90%. 


Of the 14 tests conducted during this series, five tests 
had sonic flow in both the secondary and primary nozzles. The other tests 
had subsonic flow in the primary with sonic flow in the secondary. Also 
shown listed in the table is the normalized primary nozzle flow rate. This 
term is plotted versus the chamber pressure ratio in Figure 34. 

Referring to Figure 34, the normalized flow rate increases 
to a maximum as the pressure ratio decreases. The maximum constant flow 
rate corresponds to sonic flow. The range of chamber pressure ratios over 
which primary sonic flow occurs is Pcs/Pcp The corresponding flow 

split range is .ufp > 19% and < 81%. It^^should be emphasized that these 
ranges hold only^for the particular hardware tested. Upon optimization 
of the primary nozzle expansion ratio and nozzle half angles the results 
could vary. 


The conclusions of this series of tests are: 

® Nominal primary flow to a high back pressure is 
feasible. 
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Figure 32. Dual Throat Cold Flow Comparison of Mode I and 
Mode II Aerodynamic Optimum Operating Regions 
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• OBJECTIVE 

(1) DEMONSTRATE LIMITS OF ACHIEVING SONIC FLOW IN PRIMARY CIRCUIT WITH VARIOUS 
BACK PRESSURE 

(2) ESTABLISH RANGE OF CHAMBER PRESSURE RATIOS WHICH WILL RESULT IN SONIC FLOW 
IN PRIMARY THROAT 

(3) PARALLEL FLOW PERFORMANCE 

• APPROACH 

(1) ONE PRIMARY NOZZLE POSITION - L^/R^p = 5.045 (MOST LIMITING CASE - HIGHEST BACK PRESSURE) 

(2) VARY SECONDARY CHAMBER PRESSURE FROM 10 TO 41 N/CM^ (15 TO 60 PSIA) AND PRIMARY CHAMBER 
PRESSURE FROM 38 TO 79 N/CM2 (55 TO 115 PSIA) 

(3) MEASURE SECONDARY WALL PRESSURES 




Figure 33. Parallel Flow Test Program 




TABLE VII 


PARALLEL FLOW SUMMARY 


TEST NO. 

p 

p 

CS 


Pd 

% w‘ 

s 

Wp mJjQp/P^p 

SONIC FLOW 

130 

56. 

12 

0.311 

100. 

0 

0,779 


131 

64. 

58 

0.912 

14. 

86 

0.670 


132 

66. 

58 

0.882 

14. 

86 

0.711 

Q Q' 

133 

70. 

60 

0.854 

15. 

85. 

0.747 

V o 

134 

73. 

58 

0.800 

16 

84 

0,769 


135 

109. 

46 

0.424 

29. 

71. 

0.778 


136 

65. 

58 

0.896 

13. 

87. 

0.702 

* || 

137 

64. 

57 

0.897 

12. 

88. 

0.679 


138 

81. 

56 

0.693 

18. 

82. 

0.778 

* 

139 

73. 

59 

0.810 

16. 

CO 

0.771 


140 

74. 

58 

0.788 

16. 

84. 

0.774 


141 

68. 

59 

0.862 

17.. 

83. 

0.745 


142 

81. 

55 

0.681 

19. 

81. 

0.780 

* 

143 

113. 

46 

0.407 

29. 

71. 

0.782 

* 
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III, D, Data Analysis (cont.) 


° Primary sonic flow was achieved over a range of: 

® Chamber pressure ratio <.7 
° Flow split: > 19% and < 81%. 

E. CONCLUSIONS 

' The most important conclusion from these data is that the dual 
throat concept is aerodynamical ly feasible. The performance ‘loss obtained 
with the testing hardware is as .low as 0.5 percent. The loss is minimized 
by an optimum nozzle spacing corresponding to an Af/AtS ratio of about 1.5, 
or an Le/Rtp of 3.0 for the dual throat hardware tested, while requiring 
only 4% bleed flow. The cold flow results indicate that the Mode I and 
Mode II geometry requirements are compatible and pose no significant 
design problems. 
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IV. AERODYNAMIC MIXING MODEL 


A. OBJECTIVES AND GUIDELINES 

The primary objective of this task was to develop an analytical 
model capable of predicting the aerodynamics unique with the -Mode II opera- 
tion of the dual throat thruster. 

The plume which is generated by the expanding exhaust gases 
exiting the primary nozzle and the resultant reattachment to the secondary 
wall can best be described by two processes: (1) the base pressurization 

effects caused by the reattachment of the plume to the secondary wall, and 
(2) the effects of the secondary or bleed flow (mass addition) to the super- 
sonic flow and its effect on the base pressure and plume shape. 

A survey of the literature on the above topics revealed the 
existence of a vast amount of work associated with the base pressurization 
of the wake associated with high velocity projectiles. Many models have 
been developed to predict these base pressurization effects, ranging from 
sample empirical models to highly sophisticated, finite element solutions. 
Because of the anticipated use of this aerodynamic model for repeated 
parametric studies in the future, the base pressure theories of Korst, 
Chapman and Bauer were used in the development of this model. The 
resulting model met both the resource and technical requirements of this 
aerodynamic analysis. A complete description of the model and a compari- 
son with the cold flow data are presented in this section. 

B. MODELING APPROACH AND DEVELOPMENT 

The Mode II (sustainer mode) operation is illustrated-in 
Figure 35. The primary thruster expands the flow to a supersonic condi- 
tion at the nozzle exit. Further expansion then occurs in the form of a 
Prandtl -Meyer expansion fan at the nozzle lip. The result is an exhaust 
plume with a constant pressure boundary. Flow is introduced into the 
secondary chamber so as to control the location of the plume boundary in 
order to minimize reattachment shocks which occur when the flow impinges on 
the secondary nozzle wall. The plume boundary is an exhaust streamline path 
and acts very much as a nozzle wall'. However, a shear layer develops along 
this boundary due to viscous interaction of the exhaust jet and the gases 
recirculating in the secondary chamber. Analysis of this shear layer mixing 
region and the phenomena described above are an important feature of the 
aerodynamic model. 

The assumptions employed in the aerodynamic model are listed 

below: 

° The flow leaving the primary nozzle is one dimensional and 
supersonic. 
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Figure 35. Mode II, Sustainer Mode 
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IV, B, Modeling Approach and Development (cont.) 


The exhaust is represented by an ideal gas with a constant 
expansion coefficient, y 

The shear layer is represented by the shear layer model 
developed by Chapman (Ref. 1), Korst (Ref. 2), Bauer 
(Ref. 3), and others. Assumptions implicit in this model 
have been retained. 


The shear layer is isoenergetic. 

The shear layer is treated as two dimensional.* 


The shear layer is a constant pressure surface, both in 
cross section and along its entire length. At the inner 
edge of the shear layer, the gas velocity is negligible. 

At the outer edge of the shear layer, the velocity is con- 
stant and equal to the velocity of the plume boundary 
streamline. 


The plume boundary can be represented by the methods developed 
by Herron (Ref. 4). 


Details of the aerodynamic model are discussed below. 


1 . Nozzle Wall and Plume Geometry 

Both the primary and the secondary nozzles are assumed to have 
geometries that can be described analytically by the following parameters: 

° Chamber contraction ratio. 

° Circular wall radius connecting the chamber and the 
nozzle inlet. 

° Conical nozzle inlet. 

® Circular wall radius upstream of the nozzle throat. 

° Circular wall radius downstream of the nozzle throat. 


*This assumption has been shown by Bauer (Ref. 3) to be accurate provided that 
the projected thickness of the mixing zone on the radius of symmetry is 
less than .3 of the radius of symmetry. This condition is always satisfied 
in the dual throat engine. 
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IV, B, Modeling Approach and Development (cont.) 


° Either a conical exhaust nozzle or a contour shape 
that can be represented by a spline fit. 

An exact definition of these parameters is given in the WALL 
subroutine. The inplit data write-up for the program provides the details 
concerning their use. The above parameters are sufficient to accurately 
define the geometry of a rocket nozzle. The parameters have been chosen so 
as to be identical to those used by the TDK program (Ref. 5), which is to 
be used elsewhere in the performance prediction procedure. The coordinate 
system of the nozzle geometry is centered at the nozzle throat plane, and 
the nozzle throat height is assumed to be unity. In this way, a nozzle 
can be scaled by use of only one number, i.e., the actual nozzle throat 
height (or throat radius). 

The calculations carried out by the aerodynamic model 
computer program, however, are in the "plume coordinate system". This 
system is centered at the exit plane of the primary nozzle, and the nozzle 
exit plane height is assumed to be unity. Thus, the primary and secondary 
nozzle geometries roust be converted to the plume coordinate system. This 
conversion is done automatically by the computer program. 

Once the aerodynamic model calculations have been completed, 
it is necessary to output an effective nozzle contour for the sustainer 
nozzle during Mode II operation. This is done for a nozzle operating with 
sufficient bleed flow that the plume (i.e., shear layer) boundary attaches 
just downstream of the secondary nozzle throat. Attachment further down- 
stream (or perhaps no attachment) is regarded as requiring an excessive 
amount of bleed flow. Attachment upstream of the secondary nozzle throat 
is regarded as unacceptable because of induced shock structure and aerodynamic 
losses that lower engine performance. The program thus determines a bleed 
flow rate which is a "best solution" for the operating conditions of the 
engine. 


2. Plume Boundary Shape 

A correct shape for a plume boundary caused by a supersonic 
jet expanding into a quiescent media can be obtained using the method of 
characteristics (MOC), Computer solutions of this type have been verified 
to predict plume structure for a wide range of experimental conditions. 

These MOC solutions require computer calculations that are of the order of 
several minutes per solution. Thus, they are regarded as not suitable for 
parametric analysis. Consequently, a number of simpler, correlative tech- 
niques have been developed for rapidly predicting plume shape. The basic 
approach common to these methods is to predict plume shape by matching 
MOC solutions that have either been published or calculated for the purpose 


61 



IV, B, Modeling Approach and Development (cont.) 


of providing data to match. A method of this type was selected for incor- 
poration into the aerodynamic model. The method selected to determine the 
plume boundary shape was that of Herron (Ref. 4). Before a satisfactory 
choice of method could be arrived at, a number of methods were investigated, 
programmed, and tried out in the aerodynamic model. For various reasons 
they were found to be unsatisfactory. Because of the effort expended in 
determining a proper method, the methods that were investigated are dis- 
cussed below. 


The methods investigated were those of Charwat (Ref. 6), 
Latvala (Ref. 7), Andes (Ref. 8), and Peterson (Ref. 9). The Charwat 
method was used in the original Brown Engineering Computer Program, and 
was therefore immediately available. It was found, however, that the method 
failed when applied to the operating conditions of the ALRC cold flow test 
series (Ref, 10). This problem was confirmed by an author of the OANNAF 
Plume Handbook (Ref. 11), The Charwat method was investigated for inclusion 
in the handbook, but was omitted when its limitations were found. Instead, 
the method of Latvala was recommended. 

Despite this recommendation, the method of Latvala was found 
to be highly inaccurate for a gas with an expansion coefficient differing 
from that for air (y = 1.4). For example, the plume shape for an expansion 
at Y = 1*15 presented by the Latvala method gives a plume shape much too 
expansive. The simple correction for values of y other than 1.4 that is given 
in the OANNAF Plume Handbook was found to be ineffective. This problem had 
been noted previously by Andes (Ref. 8) who developed a correlation formula 
to correct the Latvala plume shape for a broad range of conditions^ The 
Andes method, however, was found to be not applicable to the Latvala method 
as given in the OANNAF Plume Handbook and required use of a computer pro- 
gram developed by Andes. This program was found to be no longer available. 

The method of Peterson was also investigated. This method, 
which appears to be quite accurate, is contained in an appendix to 
Reference 9. It has the shortcoming of being available only for flows 
with Y = 1.15. Working out the charts on which the method is based for 
other expansion coefficients would be time consuming and would require a 
set of MOC calculations to be carried out. Because the method is essentially 
geometrical and also requires a set of charts, it is not readily suitable for 
use as a computer subroutine. Unlike the Latvala method, the Peterson 
methods use a physically correct choice of correlation parameters. 

It was thus necessary to examine other methods for pre- 
dicting plume shape. Next, the Herron method was tried and found to be 
satisfactory when suitably modified. This method uses the quantity Mb/y 
as a similarity parameter to simplify the problem. As defined in the 
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Nomenclature, Mb is the Mach number along the plume boundary and y is 
the ratio of heat capacities. Given the plume expansion angle at the nozzle 
lip, defined as the nozzle lip angle plus the Prandtl-Meyer expansion angle, 
i.e., 6 = 0 p + Av, and given the parameter Mb/y, three charts can be used to 
locate three fixed points along the plume boundary. These points corres- 
pond to x/re = 5, 10 and 18. For example, given x/re = 5, and values of 6 
and Mb/y, the charts can be used to determine a value for r/re. Herron has 
correlated this method against numerous plume shapes, ranging from 6's of 
45“ through 65° and Mb/y values of 5, 7, 10, and 13. 

Unfortunately, the available charts are so small as to be 
almost unreadable. Also, for the dual throat nozzle, the plume shapes of 
interest are for values of x/re are generally less than five, values of 
Mb/y less than five, and values of 6 less than 45°. Like the other methods 
investigated, the Herron method is intended to approximate plume shapes where 
the plume is highly under expanded and, thus, relatively large in extent. 

It was necessary, therefore, to extend the Herron method to 
the region of interest for the dual throat nozzle. This was done by 
carrying out a series of MOC plume calculations and developing a set of 
tables from the results. The plume calculations were carried out so that 
tables could be constructed over the range: 

M|j/y = 2, 3, 5, 7, 10, and 13 
6 = 30, 45, 60, 75° 
x/rg = 3, 5 

Next, it was necessary to determine a plume shape from the points obtained 
from the correlation method. It has been found that a plume boundary can 
be approximated fay a circular arc with reasonable accuracy (Refs. 7 and 8). 
Three conditions are required to determine a circle, and thus, several choices 
are possible. For example, the Latvala method, as presented in the Plume 
Handbook (Ref. 11) uses the nozzle exit coordinates, the angle 6, and the 
radius of the circle as the three conditions. A better choice in matching 
plume chape, however, was found to be three points as follows: 

(1,0), (r/rgl3,3), and (r/r^j^.B). 

In this way the plume is forced to pass through three points in the region 
of interest, although in general the plume expansion angle at the nozzle 
lip will not be satisfied. 

3. Shear Layer Model and Its Use in the Aerodynamic Model 

The shear layer model is essentially as described by 
Bauer in Reference 3. The shear layer is assumed to begin at the exit 
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of the primary nozzle where its thickness is negligible. It then develops 
as a turbulent shear layer along the boundary of the primary nozzle exhaust 
plume. Viscous interaction of the exhaust plume with gases recirculating 
in the chamber of the secondary nozzle form the shear layer. In steady 
state operation, the system will maintain a stable equilibrium base pres- 
sure. This base pressure determines the plume boundary. The theory assumes 
that a streamTtne divides the shear layer such that gases on one side of the 
streamline are recirculated into the base region. This is the so called 
"dividing streamline" concept. At steady state operation, all of the' 
gases exhausting from the primary nozzle will also exhaust from the 
secondary nozzle. Gases in the recirculating region (in this case the 
chamber of the secondary nozzle) have been trapped because they are 
degraded in total pressure to the extent that they cannot penetrate the 
static pressure rise produced by turning of the flow at the point where the 
shear layer impinges the secondary nozzle wall. It is assumed in this 
model, that the dividing streamline becomes a stagnation point at the 
secondary wall. This stagnation occurs abruptly either by a compression 
that is essentially isentropic, or by a non-isentropic compression through 
a shock structure. Details of this process govern the value that will be 
obtained for a base pressure, i.e., the pressure along the plume boundary. 
Methods by which this base pressure value can be estimated are discussed 
in Section b, below. 

Assuming that the base pressure is known from experimental 
data for the case of no bleed flow, it is possible to completely deter- 
mine the flow conditions in the nozzle, both with and without bleed flow, 
using the procedures outlined in the following steps: 

(1) Conditions at the exit of the primary nozzle are 
determined using the one-dimensional relation between Mach number and 
area ratio, i.e.. 


(2) The exhaust flow is turned at the nozzle lip and 
expands to equal the boundary Mach number M|j'. The angle of this turn, 
Av, is determined, using the Prandtl -Meyer relations for supersonic 
expansion at a point, as follows: 

Av = Vfa - Vg 
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IV, B, Modeling Approach and Development (cont.) 


where: 

v+l 1/2 9 1/2 

^ arctan -V)] - arctan [Mj^ - 1] 

1/2 1/2 1/21 
arctan (M^ - 1)] - arctan [M^^-l] 

(3) The plume boundary is determined using the modified 
method of Herron as previously discussed. Since the plume boundary is a 
circular arc, all of its geometric properties (such as its position, arc 
length to any point, and slope at any point) are known analytically. 

(4) The intersection of the plume with the secondary 
nozzle wall is determined geometrically. 

(5) The Korst method allows the shear layer properties 

to be determined exactly at the point of intersection of the plume boundary 
and the wall. All of these properties can be evaluated from various 
integrals of the error function, as described in the next section. 

(6) The plume boundary is a stream line. This stream- 
line will be turned abruptly (i.e., deflected) to follow the secondary 
nozzle wall, resulting in a shock structure. The cases of most interest 

correspond to solutions of the shock relations of the weak family. Cases 

where a lamda shock structure develops, so that the flow is shocked sub- 
sonic in part of the region, are of less interest since they correspond 
to aerodynamically poor designs. The oblique shock relations are used 

to determine the rise in static pressure across the shock and the loss of 
total pressure across the shock. Since the development of shock structure 
in the nozzle is fundamentally a two-dimensional phenomenon, no attempt 
is made to treat the flow in the nozzle downstream of the origin of the 
shock. 


(7) Once a solution has been obtained for the zero 
bleed flow case, the calculations are repeated using successively larger 
values for the base pressure. The bleed flow rate that corresponds to 
this pressure increase is then evaluated using the method presented by 
Bauer (Ref. 3) and given in the following section. The calculations 
terminate when the plume boundary is found to attach downstream of the 
secondary nozzle throat.* 


*A more exact description of steps (1) through (7), above, is presented in 
Appendix A, the problem statement for the aerodynamic model. 
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The steps described above are sufficient to determine the 
shape of the "wall" that bounds the exhaust from the primary nozzle during 
Mode II operation. In order to determine the aerodynamic performance of 
one of these solutions, it is necessary to input this wall to the TDK 
computer program. The TDK program gives a MOC solution for the nozzle. 
Shocks are ignored by TDK, but an optimum design will have weak shocks 
that can be ignored, i.e., treated as isentropic compression. The inter- 
face between the shear layer model and TDK has been automated by developing 
a spline curve fit to the primary exhaust "wall" shape, i.e., the primary 
nozzle, plume boundary, and the secondary nozzle contour downstream from 
the plume attachment point. 

a. Governing Equations for the Shear Layer Model 

A derivation of the shear layer model has been given 
in numerous reports such as References 2, 3, and 12. Much material and 
an exhaustive list of references on the subject are presented in 
Reference 13. Therefore, only a summary of the equations used will be 
given here. 


At a distance, £.j, along the shear layer as measured 
from its origin, a complete solution is available for properties transverse 
to the shear layer in terms of a non-dimensional distance, n- 


as follows: 


A dimensional distance, Y, can be obtained from n 


Y = 

where: 


CT is the "shape factor" of the shear layer. 

Both n and Y are measured from the centerline of the shear layer. The 
only value of n that is actually required by the model is no> non- 
dimensional location of the dividing streamline and this value is tabulated 
as a function of Crocco number in Reference 3. 

The velocity profile, (|>, is given by the relation 

<!> = [1 + erf o]/2. 

The Crocco number, Ca«, is also required by the method and is simply 

Ca^ = 1 - T. /T„ 

« b o 
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where: 

v-1 2 

= [1 + Y 3 • ■ 

Two integrals must be evaluated using the above quantities, and these- are 
defined as: 


(ij ^ / 


dn- 


2 2 
1 - Ca Y 


di) ^ / 


<j) dn 


2 2 
1 - Ca" <!>" 


The subscript is used to designate the upper limit of the integral. These 
integral are evaluated using Simpson's rule as the method of integration 
and a standard subroutine for evaluating the error function. Becaus.e of 
the exponential character of the error function, it is desirable to divide 
the integral ,ipto several intervals. The limit values +-«> can be 
accurately taken as + 3, which is the convention used by Bauer (Ref. 3). 


Bauer gives an expression for the mass flow between 
two streamlines that can be applied to relate the pressure and length of 
the shear layer to a bleed flow rate. 


For the case of no bleed flow, the ratio of mass 
flow entrained on the base flow side of the dividing streamline to the 
total mass flow is bound to be*: 


(G/G^) 



/Cc7(Il>nD 


If this expression is evaluated for solutions found with a larger base 
pressure, the difference in mass flow must be due to the bleed flow rate. 
Therefore, the percent bleed flow associated with these solutions is: 


*This expression differs from that given by Bauer (Ref. 3, Equation 13) in 
that the flow is assumed to be two-dimensional rather than axially symmetric, 
resulting in a somewhat simplified result. 
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where: 

i = 1 represents zero bleed flow 

i > 1 represents solutions obtained by raising the base pressure 
W 3 is the weight flow rate in the secondary chamber 

Wp is the weight flow rate in the primary chamber 

It is also necessary to use an expression taken 
from the literature for the shape factor, a. Bauer gives two such expres- 
sions. The first, due to Korst, is an empirical expression based on 
experimental data in the Mach number range from 1 to 2, which is 

a = 12 + 2.758 (Korst) 

A more basic method presented by Abramovich is pre- 
ferred by Bauer. The expression is 

o = 24 (I-]) 3 /(Iq )3 (Abramovich) 

Both expressions give the classical result that for incompressible flow 
(Mb = 0 ) the shape factor is: 


It was found that the aerodynamic model successfully correlated the ALRC 
cold flow data when the Abramovich model was used, but bleed flow vs base 
pressure was over predicted somewhat when the Korst model was used. This 
trend is because the Abramovich formulation gives lower values for a than 
the Korst method (see Figure 3 of Ref. 3).* 

Reference 13 gives data and models for about a dozen 
different relations for o. At Mach numbers greater than one, all of these 
methods give a result larger than the Abramovich model, and thus, poorer 
agreement with the data. It is somewhat surprising that the Abramovich 
model predicts the data as well as it does. The results presented in 
Ref. 13, however, were obtained in studying the base flow behind a projec- 


*Since a appears in the denominator of the expression for G/Gj, a larger 
a means that a- larger bleed flow is required to produce a given change in 
plume shape. 
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tile. Thus, the conditions existing at the dividing streamline stagnation 
point are different than found in the dual throat nozzle, or in the injector 
systems studied by Korst and Bauer. A basic difference is that the stag- 
nation point occurs in the axis for the flow behind a projectile, so that 
axi symmetric effects are critical. 

A theoretical determination of conditions that exist 
in the vicinity of the shear layer attachment stagnation point is one of 
the remaining unsolved problems in fluid mechanics. Consequently, only 
semi -empirical methods exist for predicting the pressure of the shear 
layer. This problem is discussed in the following subsection. 

b. Prediction of the Zero Bleed Flow Base Pressure 

Three options are available for determining the value 
to be used for base pressure, Pb/Po» "for zero bleed flow. 

(1) The value can be input directly to the program. 

If a measured value is available, as for example the values measured in 
Ref. 10, it should be used. 

(2) The value can be calculated automatically by 
the program. For this option, a pressure is calculated based on the ratio 
of the area of the secondary nozzle throat to the area of the primary 
nozzle throat. This value corresponds to the pressure that would be 
obtained if the nozzle were perfect, i.e., parallel flow downstream- of the 
plume attachment. 


(3) A method is also available that calculates pres- 
sure ratio by use of a correlation relating the wall angle at plume attach- 
ment to a deviation from the area ratio pressure of Method (2), above. 

The correlation is based on the cold flow results obtained during this 
program (Ref. 10). 


A more general method of determining the base pres- 
sure is desirable. This problem is discussed in detail by Roberts, ;»Ref. 12, 
where the literature on the subject is reviewed and a method is recommended. 

c. Calculation of Displacement and Momentum Thickness 
Across the Shear Layer 

The basic definitions for momentum thickness, 0£, and 
displacement thickness are 
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t 



where the integrals are taken from the dividing streamline to the interior 
edge of the shear layer. Since 


dy 



and* 


U/Ug = 


p/P£ 



1 - Ca'^cj) 

CD^ 


it follows that 


and 



r 

"D 


(1 - Cal) n - i>) 
1 - Ca^ r 

CO 



6 * = 


r n 


(1 - ca^) ^ 

1 - Ca^ <!> 

CO 


■] — dn 
a 


These expressions can be used to obtain the momentum deficit and hence, the 
thrust loss accumulated in the portion of the shear layer that exits from 
the nozzle. Note that this calculation also includes the effects of 
the bleed flow addition since this is accounted for in the shear layer 
derivation. 


*Refer to Appendix A of Ref. 12 for a derivation. 
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C. MODEL VERIFICATION 

A series of cold flow tests have been conducted by ALRC using 
dual throat nozzle hardware supplied by the Marshall Space Flight Center. 
Results of these tests were described earlier in this report. 

The aerodynamic model has been used to analyze the nozzle con- 
figurations and test conditions included in the cold flow investigation. 
Hardware configurations for the tests are shown in Figure 36. Test numbers 
given in Figure 36 are for the test series conducted during Task I of "this 
program. The hardware configuration for the previous IR&D test series was 
identical with that shown for test 101 A- 110 of the Task I test series. 

Results of each test were presented in Section III. Each 
set of tests are for a fixed hardware configuration, but with the bleed 
flow varied. Table VIII gives the separation distance, L/Rtp, between the 
exit plane of the primary nozzle and the throat, plane of the secondary 
nozzle. Also given is the bleed flow rate for each test, and the measured 
value obtained for the base pressure ratio. This information, together 
with a definition of the geometry of each nozzle wall, is sufficient to 
run the aerodynamic model computer program. 

The procedure used in carrying out the analysis consisted of 
the following steps: 

(1) The cold flow tests #727-107 through -110 were selected 
as test cases. The plotted data for these tests are in Figure 37. ' 

(2) The measured value for base pressure ratio was input to 
the program. The plume model was then calibrated so that plume attachment 
was predicted as observed. This required a multiplier of 1.8 for the 
parameter Mb/y* 

(3) Calculations were made with bothuthe Abramovich shape 
factor and the Korst shape factor (see Section 3. a. for a discussion of 
the shape factor, a). It was found that the Abramovich model was in 
closest agreement with the cold flow test results, i.e. in predicting the 
variation of Cls/(Wp + Wg) with Pb/Po- Results for two sets of cold flow 
tests are shown in^Figure 38. 

(4) Calculations were made for the remaining 34 cold flow 
tests. The plume model was not recalibrated for these calculations. 

Results obtained from the analysis are presented in Table IX, 
which shows the model predicted percent bleed flow required for plume 
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PRIMARY NOZZLE SECONDARY NOZZLE 


TEST NUMBER 

^E 

£p 

r\»T» 

's 

^s 

loiA - no 

2.244 (.8835) 

1 . 666 ^^^ 

4.562 (1.796) 

2 , 979 ( 1 ) 

111 - 116 

2.244 (.8835) 

1 . 666 ^^^ 

6.175 (2.431) 

1.626^^^ 

117 - 124 

3.092 (1.218) 

3.164^^^ 

6.175 (2.431) 

1.626^^^ 

125 - 130 

3.092 (1.218) 

3.164^^^ 

4.562 (1.796) 

2.979^^^ 


Figure 36. MSFC Contractual Test Series, Dual Throat Cold Flow 
Mode II Hardware Configurations 


TABLE VIII. - COLD FLOW TEST DATA, EACH SET OF TESTS UTILIZE A FIXED 

HARDWARE CONFIGURATION 
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107,108,110,109 

.0171 
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FIG. 19 
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.0229 
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FIG. 22 

NAS 8-32666 
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FIG. 23 
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FIG. 24 
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FIG. 25 
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Figure 38. Pressure vs Bleed Flow, Measured and Predicted Values 



TABLE IX. - PERCENT BLEED FLOW REQUIRED FOR PLUME ATTACHMENT 
AT THE SECONDARY NOZZLE THROAT 


TEST CONDITIONS 
OF FIGURE NO.: 

TEST 

PROGRAM 

TEST NUMBER 

% BLEED 
FLOW 
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FIG. 16 

IR&D 

N/A 

N/A 

N/A 

FIG. 17 
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FIG. 18 
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FIG. 23 
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FIG. 24 
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FIG. 25 
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IV, C, Model Verification (cont.) 


attachment at the nozzle throat. The agreement between this prediction 
and the cold flow measurements is quite good, which shows that the pre- 
dicted plume shape varies correctly with the bleed flow rate. Thus 
a nozzle contour, consisting of the primary nozzle wall, plume shape, and 
the secondary nozzle wall, is predicted with accuracy. 

A MOC calculation was carried out for the design condition pre- 
dicted for the test case nozzle configuration, i.e. with geometry and test 
conditions as shown in Figure 37 and with bleed flow sufficient to attach 
the plume just downstream of the secondary nozzle throat. The effective 
nozzle geometry for Mode II operation consists of the following: 

(1) the primary nozzle wall 

(2) The plume boundary predicted by the aerodynamic model 
such that attachment occurs just downstream of the 
secondary nozzle throat 

(3) the secondary nozzle wall from plume attachment to the 
exit. 

Figure 39 shows a characteristic net calculation for this geometry. 
The Prandtl-Meyer expansion fan at the primary nozzle exit can be seen in 
this figure, followed by a region of compression induced by the plume sur- 
face. This so called "plume shock" can be seen to form along a right running 
characteristic surface. At the plume attachment point a shock of greater 
strength is seen to form due to compression from the secondary nozzle wall. 

This shock is also formed by the coalescence of right running characteristic 
surfaces, and then merges with the "plume shock". A region of compression 
results that is bounded above by the secondary nozzle wall, and below 
by a right running characteristic. It should be noted that this right 
running characteristic has a positive slope that is greater than that of 
the secondary nozzle wall. Thus, most of the flow in the nozzle is 
unaffected by the shock structure. Shocks are treated by the program as 
isentropic compression. 

Figure 40 shows constant Mach number profiles for the calcula- 
tion presented as Figure 39. Since the flow calculation is isentropic, 
these are also surfaces of constant pressure, temperature, and density. 

It can be seen that the plume surface compresses the flow so that the 
Mach number is nearly constant along the plume. The divergence efficiency 
calculated at the exit of the secondary nozzle for this case was found to 
be 0.981. 


Use of these results with the performance methodology developed 
for Mode II operation is illustrated with a sample calculation presented in 
Section V of this report. 
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Figure 39. MOC Calculation for the Test Case 
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Figure 40. Constant Mach Number Profiles for the Test Case 









V. 


PERFORMANCE METHODOLOGY 


A. OBJECTIVES AND GUIDELINES 

The objectives of this task were to: (1) formulate a simplified 

JANNAF performance methodology for the dual throat thruster, and (2) utilize 
the methodology to predict the hot fire performance of a thruster of similar 
geometry to the MSFC design. 

A simplified performance prediction procedure for liquid propel- 
lant thrust chambers has been developed under the sponsorship of the JANNAF 
Performance Standardization Working Group (Ref. 16). Some of the assumptions 
in the procedure are specific to expansion nozzles of conventional design. 

In order to be applicable to the dual throat rocket engine design, the 
JANNAF procedure requires further development and modification. In the 
following, a performance methodology is developed for dual throat thrust 
chambers. This methodology represents an extension of the JANNAF methodology 
of Reference 16. It is intended to be consistent with the JANNAF methodology 
so that performance predicted for dual throat engines can be compared with 
the performance predicted for conventional engines. The computer program 
analysis tools utilized for conventional engines are also utilized for the 
dual throat concepts with appropriate input parameters as determined from 
the aerodynamic model . 

B. DUAL THROAT THRUST CHAMBER PERFORMANCE PREDICTION PROCEDURE 

1 . Method of Approach 

The JANNAF performance prediction procedures are concerned 
with prediction of the specific impulse that will be delivered by a real 
rocket thrust chamber. The approach used is to first determine an ideal 
specific impulse, IspODEs for a given rocket thrust chamber. The predicted 
specific impulse, Isppred> is then calculated by the modeling of known 
aerodynamic and thermodynamic mechanisms that degrade ideal performance. 

Ideal performance is dependent upon the physical, chemical, and thermo- 
dynamic properties of the propellant combination and its combustion pro- 
ducts, and also upon operating conditions such as mixture ratio, chamber 
pressure, nozzle expansion ratio, and ambient pressure. It is independent 
of thrust chamber design parameters, such as nozzle geometry, size, material 
and injector configurations, etc. Thus, one dimensional equilibrium per- 
formance, IspoDE provides a convenient as well as accurate reference condi- 
tion for estimating performance efficiency. 

Following the procedures described in Reference 16 pre- 
dicted specific impulse can be related to ideal specific 
follows: 

^^Ppred " ^^ER ' ’^ED ’ ’ ^^PqDE^ " 


impulse as 
Alspg|_ (1) 
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V, B, Dual Throat Thrust Chamber Performance Prediction Procedure (cont.) 


where the loss models for energy release (ER), two dimensional aerodynamics 
(2D), and finite rate chemical kinetics (K) are calculated as efficiencies; 
while boundary layer losses (BL), such as are caused by the nozzle wall 
boundary layer, are calculated directly. 

In order to develop the above loss models it is first 
necessary to establish a control surface bounding the region in which the 
analysis is to be carried out. This surface is taken as the thrust chamber 
injector face plus the chamber walls. In References 16, 17, and 18 the 
assumptions made regarding these loss models are presented and their 
interrelationships and limitations are fully discussed for rocket thrust 
chambers of conventional design. 

The dual throat chamber concept requires modification and 
further development of the JANNAF procedure if it is to be applied for per- 
formance prediction. The dual throat chamber consists of two conventional 
thrust chambers, one of which is positioned interior to the other. The 
interior chamber is referred to as the "primary" chamber, and the exterior 
is referred to as the "secondary" chamber. The engine is designed to 
operate in two modes as illustrated in Figure 41. Mode I is a booster 
mode and itself consists of two modes. In Mode I either the secondary 
thruster is operated alone (series burn), or with the primary thruster 
(parallel burn). In the series burn mode (see Figure 41a), the engine is 
of the conventional type and can be analyzed by the standard OANNAF pro- 
cedure. The ability to handle this type engine by the procedures developed 
here allows a direct comparison with the methods given in Reference 16. 

In the paraljle burn mode (see Figure 41b), the primary 
engine is also operating. This mode of operation differs from a conven- 
tional engine in that the flow in the primary chamber is expanded super- 
sonically and then shocked subsonically and then expanded through the 
secondary chamber. A variety of flow situations exist depending on the 
relative flow rates and stagnation pressures of the two chambers. 

The Mode II (sustainer mode) operation is illustrated in 
Figure 41c. The primary thruster expands the flow to a supersonic condition 
at the nozzle exit. Further expansion then occurs in the form of a 
Prandtl -Meyer expansion fan at the nozzle lip producing an exhaust plume 
with a constant pressure boundary. Flow is introduced into the secondary 
chamber so as to control the location of the plume boundary in order to 
minimize reattachment shocks which occur when the flow impinges on the 
secondary nozzle wall. The plume boundary is an exhausti streamline path and 
acts very much as a nozzle wall. However, a shear layer develops along 
this boundary due to viscous interaction of the exhaust jet and the 
gases recirculating in the primary chamber. Analysis of this shear layer 
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V, B, Dual Throat Thrust Chamber Performance Prediction Procedure (cont.) 


mixing region and the phenomena described above is accomplished with an 
aerodynamic model developed for the dual throat engine. These results 
combined with the conventional JANNAF performance procedures are necessary 
to insure proper performance prediction capability. 

2. Performance Prediction Procedure 


The dual throat thrust chamber performance procedure is 
shown in Figure 42. The procedure and each of the loss models are dis- 
cussed in the subparagraphs below. For the most part this discussion is 
limited to procedures which differ from those of Reference 16. 

The control volume for the thrust chamber has been defined 
as described previously and it is necessary to determine conditions along 
this boundary. This is done by standard procedures as is illustrated fay 
the top three rows of Figure 42. Each component of the performance model 
equation (1), is also illustrated fay a row of Figure 42. 

a. Procedure for Calculating the Theoretical Specific 

Impulse, IspoDE 

The theoretical specific impulse is calculated using 
the ODE computer program (Ref. 19). The procedure is the same as used in 
Reference 16. Because the dual throat engine concept implies two flow 
circuits, the Ispqde calculation is also divided into two flow sections 
or streamtubes. The IspooE for each streamtube is determined independent 
of the other streamtube. The total engine IspogE is determined by a mass- 
averaging the two streamtube IspooE values. 

b. The Energy Release Loss Model, her 

The energy release model is based on simplified pro- 
cedures for estimating losses due to incomplete propellant droplet vapori- 
zation and mixture ratio mal-distribution. The generalized length propel- 
lant vaporization model (Ref. 20) is suitable for this application. Input 
to the model includes propellant properties (density, surface tension, 
heat of vaporization, etc.) and injector jet diameter. The model will 
determine the amount of propellant vaporization as a function of the chamber 
shape and length. As such, it is suitable for both a design and analysis 
tool . 


The effect of mixture ratio mal-distribution for both 
the primary and secondary designs can be approximated by establishing local 
macroscale stream tubes of different mixture ratios representative of that 
produced by the injection process. This can be done through the use of cold 
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V, B, Dual Throat Thrust Chamber Performance Prediction Procedure (cont.) 


flow data from a particular injector concept or analytically from results of 
the Liquid Injector Spray Pattern (LISP)* (Ref. 21) or previous experience. 
The stream tubes, once established, are expanded from the chamber inlet to 
the nozzle exit either one-dimensional ly by a simple mass average of the 
ODE Isp or two-dimensionally using the TDK computer program with multi- 
stream tubes. 


The rate of mixing between the primary and secondary 
flows achieved during Mode I parallel burn operation can also be approxi- 
mated using a gas/gas mixing model developed by ALRC during Contract 
NAS 3-14379 (Ref. 22). The effect of this mixing process can also be 
approximated by examining the limit cases of no mixing (i-e., two stream 
tubes) and complete mixing (one stream tube). 

c. Film Cooling Losses 

If film cooling is required, the associated per- 
formance losses are predicted with the energy release model previously 
described. The flow within the film cooling region is analyzed as a 
separate flow stream from the main flow stream. These results are then 
mass averaged to obtain the total engine energy loss for both the main 
flow and the film cooling flow. 

d. The Two Dimensional Loss Model, n2D 

The two dimensional loss can be separated into a loss 
due to upstream nozzle throat curvature and a loss due to nozzle divergence 
shape. 


The loss due to upstream throat curvature affects 
only nozzle mass flow coefficient, Cp, and, thus, mass flow and C*. For 
the flow situations shown in Figure 41, nozzle Cq can be calculated by 
the following formula , — 


C 


D 



T__ (8y - 27) 

96 2304 (1 + R) 


( 2 ) 


754y^ - 757y + 3633 
276480 (1 + R)^ 


*LISP is a computer program module of the "Distributed Energy Release" 
(DER) computer program which is used to define liquid spray mass and mix- 
ture ratio distribution produced by an array of conventional rocket 
injector elements. 
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V, B, Dual Throat Thrust Chamber Performance Prediction Procedure (cont.) 


where y is averaged to_account for flow strlations (Equation 2 is rela- 
tively Insensitive to y)- 

The loss due to nozzle divergence shape affects only 
nozzle thrust. The methods by which this loss may be determined for the 
flow situations shown in Figure 41 are discussed in the following sub- 
sections. 


(1) Two Dimensional Loss, Mode I Operation 

The simplified procedure of Reference 16 gives 
charts for determining divergence loss for conical nozzles and optimum 
nozzles of a Rao type (Ref. 24). These charts are considered sufficiently 
accurate for evaluating Mode I operation. 

' (2) Two-Dimensional Loss, Mode II Operation 

The aerodynamics of the sustainer mode of opera- 
tion, i.e.. Mode II, have been described in detail in Section IV. The 
output of the aerodynamic computer program is a table of effective nozzle 
contour points suitable for direct input into the TDK computer program as 
shown in Figure 42. 


The TDK method of characteristics calculation 
is then performed to determine the effective two-dimensional expansion 
efficiency for Mode II for expansion through the primary nozzle, along 
the plume boundary extending from the primary nozzle to the attachment 
point on the secondary nozzle and finally to the exit of the secondary 
nozzle. The output from this TDK analysis includes the required n2D 
for Mode II operation. For this case, shocks are considered as i sen- 
tropic compressions. 


In the event that the above procedure is con- 
sidered to be of insufficient accuracy, a method of characteristics cal- 
culation with an oblique shock capability is required. This more compli- 
cated treatment might be required if the nozzle is poorly designed such 
that plume impingement upstream of the throat gives relatively strong shock. 
In this case, if the secondary nozzle expansion ratio is large, the simpli- 
fied TDK calculational path might again be inadequate due to important 
shock losses. 


e. The Chemical Kinetics Loss Model, 

The chemical kinetics loss model is identical to that 
of Reference 16. The high chamber pressures planned for the dual throat 
engine make this loss negligible except possibly at high area ratio. Several 
sources of generalized kinetic loss charts and figures exist for LOX/H 2 and 
LOX/RP-1 propellants which can be used to approximate the kinetic losses. 
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V, B, Dual Throat Thrust Chamber Performance Prediction Procedure (cont.) 


A model for each of these systems has been automated in a computer sub- 
routine called CHEW (see Appendix C). An exact solution can also be found 
for the dual throat configuration through the use of either the ODK or 
TDK computer models. 

f . Boundary Layer Loss Model , ng[^ 

The nozzle wall boundary layer loss is to be calcu- 
lated using the BLPL computer program developed at the Aerojet Liquid 
Rocket Company. This computer program is an implementation of the turbulent 
boundary layer chart procedures given in Reference 17. Input to the pro- 
gram includes operational parameters such as Pc, wall temperature, and y, 
design parameters such as Rt> e» and nozzle shape (i.e., conical or bell). 
The output includes an estimated boundary layer thrust decrement and 
displacement thickness. Other more rigorous computer solutions also exist 
including the BLIMP-J and TBL Programs. In addition, the aerodynamic 
computer model calculates the momentum and displacement thickness growth 
occurring within the shear layer during Mode II operation. This informa- 
tion can be used with the BLIMP or TBL computer codes to determine the 
total effective boundary layer loss for Mode II operation. Using this 
procedure, the effect of the bleed flow on overall performance is 
accounted for with the shear layer model in terms of an overall boundary 
layer performance loss. 

If required, the transpiration cooling of the nozzle 
wall produces another loss which must be evaluated. The transpiration 
cooling mechanism which is modeled as mass injection into the boundary 
layer is a special case of the boundary layer loss model and the resultant 
loss caused by mass injection into the boundary can be predicted with the 
BLIMP-J computer program. 

C. DUAL THROAT THRUST CHAMBER PERFORMANCE PREDICTION - 

SAMPLE CASE 

A sample performance prediction case using the procedures des- 
cribed in Section V.B. is presented in this section. The thrust chamber 
geometry for this sample case was chosen to be identical to the dual throat 
cold flow test hardware. By choosing the geometry in this manner, both 
a sample performance calculation and a cold flow to hot fire comparison were 
made available. The important geometric dimensions are given in Hgure 43. 

The thrust chamber conditions and the predicted performance 
values for this sample case are presented in Table X. The one-dimensional 
equilibrium specific impulse values were obtained using the CHEM computer 
program described in Appendix C. As indicated in Table X, this performance 
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LE X. - SAMPLE CASE PERFORMANCE PREDICTION 






MODE I 

MODE 

II 


PARALLEL 

SERIES 




PRIMARY 

SECONDARY 

SECONDARY 

PRIMARY 

SECONDARY 

PROPELLANT 

LOX/LHg 

LOX/RP-1 

LOX/RP-1 

LOX/LH^ 

EXHAUST 

Pc,N/cm^ 

1,379 

689 

689 

1,379 

80 

(psia) 

(2,000) 

(1,000) 

(1 ,000) 

(2,000) 

(116) 

0/F 

7.0 

2.8 

2.8 

7, 

7. 

Tc “K 

3,680 

3,701 

3,701 

3,680 

3,680 ‘ 

(“R) 

(6,624) 

(6,662) 

(6,662) 

(6,624) 

(6,624) 

HEAT CAPACITY 
RATIO. Y 

1.18 

1.13 

1.13 

1.18 

1.18 

C*. lt)/s 
(ft/sec) 

2,246 

(7,370) 

1,782 

(5,848) 

1,782 

(5,848) 

2,246 

(7,370) 

2,246 

(7,370) 

ft kg/s 

5.830 

17.957 

28.289 

5,830 

.321 

(LBin/sec) 

12.852 

39.588 

55.752 

12.852 

.708 

% OF TOTAL 

24.5 

75.5 

100 

94.78 

5.22 

AREA RATIO, 

2.979:1 

2.979:1 

2.979:1 

20.5:1 

20.5:1 


MASS AVERAGED 


MASS AVERAGED 

IspooE.sec 

305.5 


287.0 

431.0 


nk 

t)D 

nER 

Aisp BL.sec 
Delivered Isp.sec 
% Isp ODE 
Thrust, N 
(LBf) 

.9931 

.9918 

.9950 

1.8 

297.6 

97.4 

69419. 

(15606.) 


.9913 

.9918 

.9950 

1.8 

279.0 ' 

97.2 
69192. 
(15555.) 

.9988 

.9805 

.9950 

7.14 

412.8 

95.8 

24901. 

(5598.) 

.9880* 

4.62* 

418.6* 

97.1* 

25248. 

(5675.)* 


* SIMILAR VALUES FOR CONVENTIONAL, 20:1, 12,5“ CONICAL 


V, C, Dual Throat Thrust Chamber Performance Prediction - Sample Case (cont.) 


from this nonoptimized hardware is about '\-97% for Mode I and '^ 96 % for Mode II. 
The Mode II efficiency is only 1.3% less than that possible from a geo- 
metrically similar, conventional engine which compares very well with the 
approximately 1% loss obtained with the cold flow tests at the same Le/Rtp 
of 2.883 {Test 727-110). 
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VI. PRELIMINARY ENGINE SYSTEM MODEL 


A. OBJECTIVES AND GUIDELINES 

The objective of this task was to define how the analytical 
results could be incorporated into an engine system model. Since testing 
was required on the contract to more fully determine the aerodynamic model, 
the scope of the engine system model task was greatly reduced. 

B. ENGINE SYSTEM DEFINITION 

A logic chart is given in Figure 44 showing how the results 
from this contract, combined with related study results, will lead to 
the definition of a baseline dual throat engine system. The initial step 
is to establish a preliminary baseline engine (Figure 45 and Table XI) from 
which parametric data can be generated. 

Preliminary designs and parametric data for other engine concepts 
are available based upon a reference thrust level of 2.70 MN (607,000 Ibf) 
from the Advanced High Pressure Engine Study (NAS 3-19727), and vehicle 
mission studies were conducted utilizing engines of this thrust level. 
Therefore, a dual throat engine definition analysis could take advantage 
of all this prior work if the reference point were taken as 2.70 MN. 

The typical preliminary baseline engine given in Table i XI has 
a thrust split Fi/Fg' =2.42 calculated based on the ratio of vacuum thrusts. 
This split resulted from a selection of a 60% LO 2 /RP-I and a 40% LO 2 /LH 2 
sea level (Mode I) thrust contribution. This percentage during Mode I 
operation was found to be optimum in recent studies at NASA/Langley 
Research Center (Ref. 25) utilizing the dual expander engine. It is 
expected that the dual throat engine design will optimize at a different 
point than the dual expander engine. The similarity of the two concepts, 
however, makes the selected split a good starting point for the conduction 
of parametric studies. 

C. ENGINE SYSTEM MODEL 

The major emphasis of this contract has been to generate an 
aerodynamic model and to incorporate this into a performance model for 
the dual throat thruster. To obtain an engine system model, the per- 
formance model needs to be incorporated as a subroutine with computer 
models providing cycle power balance, heat transfer/structures/life 
determination, turbomachinery design, thrust chamber design, nozzle 
design, component weight scaling, and mission sensitivity. The bulk of 
this requirement can be achieved in a three task effort as outlined in 
Figure 44 and given in more detail in Figures 46-48. 
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{ Figure 44. Dual Throat Engine Definition Logic Chart 























Figure 45. Preliminary Baseline Dual Throat Engine 



I TABLE XI. - PRELIMINARY BASELINE DUAL-THROAT ENGINE 


thrust (SL) MN (lb) 

Thrust (V) MV (lb) 

Mixture Ratio 

* Chamber Pressure N/CM (psia) 
! Area Ratio 
: ODE Is (SL) S 
; ODE Is (V) S 
: Ig Efficiency % 

I Is (SL) Deliv S 
Is (V) Deliv S 
Total Flow Rate Kg/S (Ib/S) 
Fuel Flow Rate Kg/S (Ib/S) 
Oxidizer Flow Rate Kg/S Ib/S 
C* M/S (ft/S) 

Throat Area CM^ (in. 2) 

’ Exit Diameter CM (in.) 



F^/Fg = 2.44 


60% X 1 
LO 2 /RP-I 

1.62 (364,200) 
1.90 (427,400) 
2.8 

1724 (2500) 
(50:1) 

310.7 

364.6 
97 

301.4 

353.7 

548.1 (1208) 
44.2 (318) 

403.9 (890) 
1803 (5914) 

573 (88.9) 

191.0 (75.2) 


40% X 1 
LO2/LH2 

1.08 (242,800) 
1.23 (277,300) 

7.0 

2068 (3000) 
(50:1) 

395.9 

452.2 
98 

388.0 

443.2 

283.9 (626) 
35.5 (78) 

248.4 (548) 
2255 (7399) 

309 (48.0) 

140.4 (55.3) 


Mode 1 

LO2/RP-I & LH2 

2.70 (607,000) 
3.13 (704,700) 
3.63 

1724/2068 

50:1 


97/98 

330.9 

384.2 

832.0 (1834) 

179.7 (396) 

652.3 (1438) 

883 (137) 

237.1 (93.3) 


Mode 2 
LO2/LH2 


1.28 (288,700) 
7.0 

2068 (3000) 
142.6:1 

470.8 
98 

461 .4 

283.9 (626) 
35.5 (78) 

248.4 (548) 
2255 (7399) 

309 (48.0) 

237.1 (93.3) 
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Figure 46, Task I - System Evaluation 
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Figure 47. Task II - Parametric Data 
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VII. CONCLUSIONS AND RECOMMENDATIONS 
A. CONCLUSIONS 

The dual throat cold flow test program has been successfully 
accomplished with the completion of 69 tests on IR&D and contract funding. 
The analysis of these data has resulted in a better understanding of the 
performance trends and characteristics of the dual throat engine. 

The cold flow test results coupled with the aerodynamic base 
flow theories of Korst, Chapman and Bauer form the foundation of the 
aerodynamic model developed in this program. Verification of the model 
with the cold flow data gives validity to the analytical techniques 
used. 


The performance prediction methodology developed for the 
dual throat engine is based on the JANNAF prediction procedures. The 
aerodynamic model provides the only additional information needed in 
conjunction with the traditional performance codes and logic given in 
the JANNAF procedures. Because this methodology is consistent with 
the JANNAF procedures, a one-to-one comparison between the dual throat 
and conventional engine performance is possible. Comparisons of this 
type will prove to be invaluable when determining the advantages and 
disadvantages of the dual throat design. 

B. RECOMMENDATIONS 

The aerodynamic model should be modified to predict the base 
flow pressure using the iterative method given by Roberts, Reference 12. 

It will probably be necessary to modify the method so that agreement can 
be obtained with the cold flow data. This method predicts base flow 
pressure directly as a function of bleed flow, in contrast to the indirect 
method currently used. Either a direct or indirect method can be used once 
the zero bleed flow solution has been obtained. 

The isoenergetic shear layer model should be modified to treat 
the mixing of non-isoenergetic streams; i.e., the case where the bleed 
flow has a different stagnation temperature than the primary flow. The 
situation where the bleed flow has a different expansion coefficient than 
the primary flow should also be treated. 

Plotting should be added to the program. Plots of interest 

are: 


(1) wall -plume-wall contours 

(2) base pressure vs bleed flow 

(3) shear layer profiles 
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VII, B, Recommendations (cont.) 


A method should be developed for obtaining two dimensional 
axially syiwnetric flow field solutions that include the effect of shock 
waves. The problem can be important when; 

(1) The operating conditions are off-design so that shock 
structures develop. It is important to know the effect 
this can have on performance. 

(2) Mode II operation at design condition will necessarily 
contain shock waves. Although negligible at low expan- 
sion ratio, the shock(s) can become strong downstream. 

A method should, therefore, be developed for estimating 
the effect of non-isentropic shock losses on performance. 

It is also important to be able to predict the nozzle wall 
pressure so that the possibility of flow separation can be 
avoided. 

Since shock wave development in a nozzle is fundamentally a 
two-dimensional problem, a two-dimensional computer analysis must be used. 

It is recommended that the NAP Code (Ref. 14) be applied to one or more nozzle 
configurations to determine if it is cost effective for this type of prob- 
lem. A second possibility is to use an existing MOC program that either 
treats, or can be modified to treat, shock structure. Unlike the NAP 
analysis, this approach will not work for flows with lambda shock structure. 

A series of calculations should be carried out using the aero- 
dynamic model and the performance prediction model that will provide design 
criteria for the dual throat nozzle optimization. A MOC program would be 
used to establish optimum nozzle contours for selected designs. Both 
Mode I and Mode II performance need to be considered. For Mode Ir only the 
secondary nozzle contour is involved and the optimum contour can be deter- 
mined by the conventional Rao method (Ref. 15). For Mode II, both the 
primary and secondary nozzle contours are involved. These contours also can 
be determined using the method of Rao and supersonic start conditions deter- 
mined -from the TDK analysis of the primary nozzle and plume. 

For the primary nozzle a suggested approach would be to calcu- 
late a family of Rao nozzle contours, each of fixed expansion ratio, but of 
variable length (i.e., the lip angle constraint would be removed). A con- 
tour would be chosen from the results by trading off lip angle vs length. 

For the secondary nozzle, a Rao contour would be calculated 
that is optimum for Mode II. A standard optimum nozzle contour would be 
generated, i.e., a nozzle that is optimum performing for its length and 
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VII, B, Recommendations (cont.) 


of fixed area ratio. However, the calculation would be started from 
supersonic conditions calculated by the TDK program using the contour of 
the primary nozzle, extended by the plume boundary. The result will be 
the design of two different contours for the secondary nozzle: (1) a 

nozzle optimized for Mode I operation, and (2) a nozzle optimized for 
Mode II operation. 

Calculating the performance of each nozzle in each mode of 
operation will allow a trade-off to be made such that an optimum contour 
can be selected for the mission requirements. 

The aerodynamic model should be correlated with hot fire test 
data. Data at actual engine conditions could be used to verify and fine- 
tune the aerodynamic model and performance methodology presented herein. 
Likewise, valuable heat transfer data would be made available during such 
tests. 


The CHEM program described in Appendix C should be upgraded 
for greater flexibility. The present tables should be extended to include 
additional propellant combinations (currently LOX/RP-1 and LOX/LH 2 are 
available). Wider ranges for Pc, area ratio and mixture ratio are also 
needed to provide adequate coverage of anticipated analyses. 
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APPENDIX A 


AERODYNAMIC MODEL FOR MODE II 
OPERATION OF THE DUAL THROAT ENGINE 


A ste'p-by-step procedure is given below that describes a computer 
program that models the aerodynamic behavior of the dual throat engine 
when operating in the sustainer mode, i.e.. Mode II. The FORTRAN computer 
program is called BASE because the aerodynamic model is constructed on a 
theory for supersonic base flow. The main subprogram for the BASE code is 
called MAIN and each step discussed below is labeled with comment cards 
as a separate paragraph in the FORTRAN code for MAIN. 
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step 1: Default input values and program constants are set- The pro- 

gram data is read, first a one card lablej, and then the NAME • 
LIST data, i.e.-, $ DATA. The above input is printed. Program 
constants are calculated from the input and input angles are 
converted from degrees to radians- The following quantities 
are then calculated: 


Step 2: 



V 


P 




) 


X* = (Le/r*)/r 

^ P p 


•If IwALL = 2, the primary and secondary nozzle walls are 
generated. Subroutine WALL is used to calculate vectors: 


X, y, and dy/dx 

for each nozzle. These vectors are generated in a coordinate 
centered at the throat of each nozzle and normalized by the 
throat radius. Next, each vector is transformed into a 
coordinate system centered at the exit of the primary nozzle 
and normalized by the primary nozzle exit plane radius. This 
is called the "plane coordinate system". 

The transformations are: 


primary nozzle; 
x_ 


(x - X )/r , where Xg is 

P P P 

routine WALL 
to y^ . 


provided by sub- 
and corresponds 




dx 
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step 3: 


Step 4: 


secondary nozzle; 


X 


s 




r*r 


P e. 


+ O-e/--;) 





Either Mgp is input directly, MEPI ^1; or if MEPI < 1, Mgp 
is found as a function of sp using subroutine MACHS, then: 


/P. 

6p 0 


[1 + ^ Mg ] “ 




[1 + ^ m2 ] 


-1 


r„. = r? 




then 


(Vo' = (Vo’, 
Pb' ■ 


otherwise 


and 


;Y-1 




input 

n 
-1 


i/a 


■=b ' 


“ ^b^%^ using subroutine 


MACHS. 
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step 5: 


Step 6: 


Y-1 2 

(p^/Po) = n * Y i-i 


i = 1 


( VT„) = [1 + ^ H^] 


-1 


Values are printed for Pjj/P-j» ^5/^0 

2 

The Crocco Number, Ca , is calculated as 

00^ 

Ca^ = 1 - T. /T„ 

CO D O 

and the following are calculated and printed: 

Up from subroutine ETADS 

(I,) from subroutines IDNES 

I rip 

(Ig) from subroutine IZERDS 

(It) from subroutine IDNES 

‘ 3 

a from subroutine SIGMAS 



subroutine PLUME is used to find v. , v , x-, y., y. , s. 

^ w J J u 

subrountine MAX is used to find the index corresponding to 
the maximum plume radial corrdinate. (i.e.. The maximum 
plume length corresponds to this coordinant.) 

The intersection of the plume xj, yj with the secondary nozzle 
wall is found. This point is called 

X, y, y , ?, and ® 

_ _ s 

where (x, y) are thejDoint coordinates, y is the slope of the 
wall at this point, s is the arc length of the plume, and 

{^) is the slope of the plume at this point. 
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step 7: 


Mass flow ratio is calculated as 


Step 8: 


(G/G^) 


(i) _ 


Y+1 

9 1 2 Ph 2(i-l ) 

p 0 


Ca‘ 




The percent increase in mass flow over zero bleed flow is 
calculated as 



0 for 1=1 

C(G/G^)^^^ - (G/G^)^^^]*100 for i>l 


Values of the intersection point properties calculated in Step 6 
and (P5/P0) are saved. Next, the pressure is incremented 


(W = (Pb/P„) + AP 

- X- L 1/2 

"b - - 1]) 

i = i + 1 

If any of the following conditions are met, go to Step 9: 
i > i 

max 

or 

(P. /P^) > (P^ /P„) 

bo ' 6p 0 ' 
or 

> WSWP(NW), which is the largest value 
requested for output (see $ DATA). 

Next, if I WALL = 1 (i.e., cylindrical wall option) go to 
Step 5. 

"k 

If x < X2 go to Step 5. 

Otherwise, proceed to Step 9. 
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step 9: 


Using Wj/Wo as the independent variable, interpolate in the table 
saved in Step 8 using as the fndependent variable = 0, 

WSWP(l) WSWP(NW) to obtain values for x, y, 



Print the following: 


\f’ ^w’ ®w’ ^b ’ ^tJ^o 


where 


W /W are the input values 

b P 


w 


w 


'w 


w 


e° 

w 




plume/wall intersection in plume 
J coordinates 


= y 


r*r 


Tx - a /r*l -i-1 


I plume/wall inter- 
) section in secondary 
wall coordinates 


^ arctan 


wall angle at plume/wall 
intersection 


180 


arctan y‘ 


plume angle at plume/wall 
intersection 


a - e. 


streamline deflection angle 
corresponding to induced shock 


Pb/Po 


pressure ratio along plume 
boundary 




pressure ratio behind the shock 




ratio of total pressure behind the 
shock to the total pressure 
before the shock 
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The angles e°, a°, and 6° are shown in the figures below; 



Thus, the plume streamline is deflected through an angle, 
65. This deflection angle corresponds to a shock angle, 
0s, as shown in the figure below: 



Subroutine WAVANG is used to determine the shock angle, 65; 
the static pressure ratio across the shock, P2/Pi» 
total pressure ratio across the shock, P-p /?-n (same as 

'2 'l 

Pt /P«» above). 

J 2 0 

The above calculations are also carried out and the results 
printed for the "best solution found". This solution corresponds to the 
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bleed flow found to be just sufficient to cause the plume to attach just 
downstream of the secondary nozzle throat; i.e., it is the last solution 
found as determined by the criteria given in Step 8. 

Next, subroutine FIT is used to construct a spline fit of the 
effective nozzle contour corresponding to the "best solution found". This 
contour consists of the primary nozzle wall, the plume, and the secondary 
nozzle wall from the plume attachment to the exit. The coordinate system 
is that of the primary nozzle. Results are output to be input directly 
into the TDK computer program. 

The momentum thickness, eg, and the displacement thickness, 6*, 
are calculated at the shear layer cross-section at plume attachment. The 
integrals -are taken from the invisid edge of the boundary layer to the 
dividing streamline; i.e., stagnation streamline, using subroutines IDELTS 
and ITHETS, respectively. The results are printed. 

The program then proceeds to Step 1 to process the next case. 
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APPENDIX B 


PROGRAM INPUT, AERODYNAMIC MODEL 


CARD 1: 


The 1st card to be input is used to provide a printed heading for 
the computer output. Columns 2 through 72 are available for text. This 
card must always precede each case. The remainder of the input data is 
read using NAMELIST with the name $DATA, as described below. 

$DATA 


Geometry; Option 1, Cylindrical Secondary Wal 1 , see Figure 49. 


INPUT ITEM 


DEFINITION 


IWALL = 1 , 

RS 

RSP 

THETAP 

EPSP 


The secondary nozzle wall is a cylinder 
r*, radius of the secondary nozzle wall 
r*p, radius of the primary nozzle 
0p, exit cone half angle, primary nozzle 
Ep, expansion ratio of the primary nozzle 


The items r* and r*p, above, can be input with units of inches, feet, cm, 
etc. The input item ep, above, is used only to determine the Mach number 
at the exit of the primary nozzle (not required if MEPI is input). 
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Figure 49. Nozzle Geometry for Option 1 



Figure 50. Nozzle Geometry for Option 2 
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Geometry; Option 2, Full Geometry Option, see Figure 50. 


IWALL = 2, 

Complete wall geometry specified for both 
primary and secondary nozzles 


Primary Nozzle 

ECP 

£(. , primary nozzle chamber contraction ratio 

RIP 

p 

Rip, wall radius connecting chamber and inlet 
or primary nozzle 

THIP 

0 I„, inlet angle for primary nozzle 

RWTUP 

Rjj , wall radius on upstream side of the 
throat, primary nozzle 

RSP 

r*p, throat radius of primary nozzle 

RWTDP 

Rdjj, wall radius on downstream side of the 
throat, primary nozzle 

THETAP 

9p, exit cone half angle, primary nozzle 

EPSP 

e , expansion ratio of the primary nozzle 

r 


Secondary Nozzle 

EC 

Ec, secondary nozzle contraction ratio 

RI 

Rj, wall radius connecting chamber and 
inlet of primary nozzle 

THETAI 

e^j, inlet angle for secondary nozzle 

RWTU 

Rm, wall radius on the upstream side of the 
throat, secondary nozzle 

RS 

r*, throat radius of secondary nozzle 

RWTD 

Rfj, wall radius on downstream side of the 
throat, secondary nozzle 

THETA 

e°, exit cone half angle, secondary nozzle 

EPS 

e, expansion ratio of the secondary nozzle 
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Geometry; Options 3, 4 or 5. Full Geometry Option with Spline Fit 
Exhaust Contours - 


For the spline fit wall contour options, contour points are input too 
replace the conical exit cone of either the primary nozzle, secondary nozzle, 
or both. Other inputs are the same as in geometry Option 2, described above, 
except that the inputs for en, and/or e are not required. The exit cone 
half angle is interpreted as^the spline contour attachment angle. The spline 
options are as follows: 


I WALL = 3, 
IWALL = 4, 
IWALL = 5, 


Primary Wall 


Secondary Wall 


Spline 
Cone 
Spl i ne 


Cone 

Spline 

Spline 


The spline coordinates are input as follows: 


X4P(2)* = 

xpj, axial coordinates for primary wall 
spline contour 

Y4P(2) = 

yp. , radial coordinates for secondary wall 
spline contour 

N4P 

np, i =1, 2, ... np above, np £ 20 

TH4P 

e°ep, spline contour exit angle for pri- 
mary nozzle, degrees 

X4(2} = 

Xi, axial coordinates for secondary 
wall spline contour 

Y4(2) = 

y-j, radial coordinates for primary wall 
spline contour 

N4 

n, 1,=1, 2...n above. N £ 20, 

TH4 

e®e» spline contour exit angle for secondary 


nozzle, degrees 
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Relative Positioning of Primary and Seconda'ry 
Nozzles (IWALL > 2 Cases Only) 


XLE ' Le/r*p, axial distance from the exit plane 

of the primary nozzle to the throat plane 
of the secondary nozzle, normalized by 


*The first point on a spline contour is automatically calculated by the 
program from the exit contour half angle, ep®, or 0°, thus each of these 
input vectors begins with item 2. 
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FLOW PROPERTIES AND OPTIONS 


INPUT ITEM 
DP 

FCTR 

G 

IPNCH 

MEPI 

NW 

OPSIG 

PBPO 

WLPRN 
WSWP (1) 


DEFINITION 


ASSUMED 

VALUE 


Pressure ratio increment for generating bleed .002 

flow v.s. pressure table. 


Factor for scaling plume shape. Used as a 1. 

multiplier for the parameter Mj^/y in the Herron 
method. 

Y, expansion coefficient 1.4 

If IPNCH = 1, geometry inputs for TDK will be 0 

punched in the NAMELIST format . 

If IPNCH = 0, no punched output, only printed 
output. 

M . (REAL type) exit Mach No., primary nozzle. 0 


If M^ <1 then values will be calculated from e 
hp P 

assuming one-dimensional flow 

n , number of entries in WSWP 13 

w 

n„l50 

Option for Calculating a 1 

If OPSIG =1, 0 from Abramovich 

If OPSIG = 2, 0 from Korst 

Pressure ratio for the plume boundary for zero 0. 

bleed flow. If not input, the program will calcu- 
late a value corresponding to the 'area of the 
secondary nozzle throat divided by the area of the 
primary nozzle throat using the relationship for 
one-dimensional flow. 


If a value other than zero is input, print of the 
secondary and primary nozzle wall tables will be 
suppressed (to be used on successive cases when 
the wall -geometry does not change.) 

(Wg/^lp)^- * 100, for a given value of P[j/Pq 

percent bleed flow output will be 
given for these values, i = 1, — , n^ 


0 . 


0,1, 2,3, 4, 5, 

6,7,8,9,10, 

12,15 


$END 
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APPENDIX C 


THE CHEM COMPUTER PROGRAM 


This program calculates equilibrium and kinetic performance values 
for Isp and C*. The propellants must be either LOX/RPl or LOX/H2. 

For LOX/RPT systems the chamber pressure must be between 300 and 1000 
psia, the mixture ratio between 2.4 and 3.4, and the nozzle expansion ratio 
between 1 and 300. Chamber enthalpies for LOX and RPl are assumed to be 
-3093 and -6200 calories/mole, respectively. 

For. L0X/H2 systems the chamber pressure must be between 300 and 
2000 psia,' the mixture ratio between 5 and 10, and the nozzle expansion 
ratio between 40 and 3000. Chamber enthalpies for LOX and H2 are assumed 
to be -3093 and- -2154 calories/mole, respectively. Kinetic value for Igp 
are calculated based on an assumed throat radius of 2.289 inches. These 
values are corrected as a function of throat radius as follows: 



KIN 


I -log (2.289/r*). 

^Pkin 

2.289 


Program Input 
$DATA 

FUEL = The propellant is either LOX/RPl (=0) or L0X/H2 (=1). 

XMAR = Propellant mixture ratio 

PC = Chamber pressure, psia 

EPS = Nozzle expansion ratio 

RT = Nozzle throat radius, inches 

SEND 


Program Output 

The following items are output: 


^^PkIN’ 

C*od£. ft/sec 


C*Kim, ft/sec 
"^KIN 
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